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FOREWORD 


Simulating  the  interaction  of  Lerrain,  vehicle,  and  man  is 
essential  Lo  improving  land  mobility  technology  —  a  goal  of  the 
U.S.  Army  Tank-Automotive  Command  (TACOM)  and  the  U.S.  Army  Corps 
of  Engineers  WaLerways  Experiment  Station  (WES).  Under  TACOM 
contract,  Grumman  has  been  actively  supporting  this  effort  and 
has  just  completed  the  development  of  a  mathematical  model  for 
pneumatic  tires-soil  interaction  an  essential  element  in  the 
terrain-vehicle-man  system  —  which  is  Lhe  subject  of  this  report. 

For  Grumman,  participation  in  this  program  began  in  1971, 
when  the  TACOM-WES  first  generation  terrai n-vehic lc-man  simulation 
called  the  "AMC  :71  Vehicle  Mobility  Model"  was  completed.  This 
model  made  it  clear  that  one  sub-model  needing  improvement  was  the 
simulation  of  a  wheeled  vehicle  moving  in  suit  soil,  and,  there¬ 
fore,  TACOM  contracted  with  Grumman  to  deveion  Li  computerized 
mi  del  of  rigid  wheel-soil  interaction  or.  the  basis  cf  the  plas¬ 
ticity  theory  for  soil.  The  drawbar  pull,  torque  requirements, 
slip,  and  sinkage  for  rigid  wheel  calculation  by  the  model  was 
validated  in  an  experimental  program. 

A  follow-on  program,  the  results  of  whit,  h  are  presented  in 
this  report,  began  in  1972.  Here  the  interaction  concept  is  ap¬ 
plied  to  pneumatic  tires,  and  a  mathematical  mode,  is  developed 
for  computing  Lire  performance  parameters  from  tjre  inputs  (di¬ 
ameter,  width,  inflation  pressure,  and  defied  i^n)  and  from  soil 
inputs  (cohesicn,  friction  angle,  and  bulk  density).  This  mathe¬ 
matical  model  is  the  first  step  toward  the  development  of  a  Lire 
performance  sub-model  in  the  AMC  Vehicle  Mobility  Model  that  can 
predict  the  ;pe_d  of  wheeled  vehicles  in  soft  soil,  accurately 
assess  torque  and  fuel  consumption  requirements  in  specific  mis¬ 
sions,  and  he  suitable  for  eventual  application  in  a  vehicle  dy¬ 
namic  mode  1  . 
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ABSTRACT 


The  role  of  tire  deformation  in  tire-soil  interaction  is  dis¬ 
cussed,  and  a  mathematical  model  is  described  that  incorporates 
the  qualitative  relationships  obtained  from  experimental  data  de¬ 
scribing  tire  deformation  as  a  function  of  the  relative  stiffness 
of  the  tire  to  the  soil.  The  mathematical  model  assumes  that  the 
centerline  geometry  cf  tires  consists  of  straight  lines  and  Icga- 
ithmic  spirals  and  that  tnere  is  a  certain  limit  to  the  normal 
pressure  that  may  develop  beneath  a  tire.  This  limiting  normal 
pressure  is  a  tire  property  and  is  dependent  on  tire  inflation 
pressure.  Soil  is  modeled  by  its  Coulomb  strength  paiameters. 
Plastic  equilibrium  conditions  in  soil  determine  the  interface 
stresses  where  the  normal  stress  is  less  than  the  limiting  normal 
pressure . 

Values  of  free  parameters  in  the  model  that  determine  tire 
centerline  geometry  were  analyzed.  Tire  performance  prcaictions 
were  compared  wi^h  experimental  data  for  a  large  number  of  com¬ 
bination  of  the  free  parameters,  and  relationships  were  established 
for  the  estimation  of  these  parameters  from  tire  deformation  mea¬ 
surements  on  a  rigid  surface.  Relationships  between  slip  and  mo¬ 
bilized  interface  shear  stresses  were  established.  On  the  basis 
of  the  tire-soil  model  a  computer  program  was  developed  that  com¬ 
putes  the  pull  and  torque  coefficients  and  the  sinkage  for  a  given 
load  and  tire  from  the  following  inputs:  slip,  inflation  pressure, 
Lire  detlection  on  rigid  surface,  tire  dimensions  and  soil  cohesion, 
friction  angle,  and  unit  weight.  Alternatively,  cone  penetration 
data  may  be  used  with  purely  frictional  or  purely  cohesive  soils. 
Predictions  by  the  tire-soil  model  are  validated  for  available  ex¬ 
perimental  da  ta  . 
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angles  defining  start  and  end  of  deflection 
central  angle  at  end  of  forward  field 
central  angle  at.  end  of  rear  field 
entry  any  Le 
rear  angle 

angle  of  tangent  of  log  spiral  with  normal  Lo 
the  radius 

coefficient  in  Eq.  (1) 

unit  weight  of  soil  (lb/cu  i>) 

angle  of  inclination  of  resultant  sLress  to 
normal,  angle  of  shear  mobilization,  deflection 
in  Waterways  Expo.,  intent  Station  terminology 

deflection  coefficient  in  Eq.  (2) 

angle  between  x  axis  and  major  principal  stress 
J  -  21  /d 

"A  -  :/2 

tolerance  limit 


Oj  +  -j )  /  2  + 

normal  stress 


principal  stresses 
shear  stress/shear  si  .  m 


maximum  available  shear  strength 
mobilized  shear  strength 
angle  of  internal  friction 
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I.  SCOPE  OF  WORK 


i' he  si'iipi'  o!  work  as  described  in  the  REP  work  statement  for 
this  contract  was  the  application  of  plasticity  theory  foi  soils 
to  the  nroblem  of  pneumaLic  Lire-soil  interaction  and  the  develop¬ 
ment  of  a  mathematical  tire-soil  model  for  tin*  purpose  of  Lire 
performance  predictions.  Within  this  general  scope  emphasis  was 
placed  on  the  following  items: 

Application  of  plasticity  theory  methods 
to  the  calculation  ol  stresses  it  a  de¬ 
formed  tire-soil  interface 

fl  Formulation  of  a  mathematical  tire-soil 
mode  l  for  the  description  of  tire-soil 
i  nte ’"action 

c  Development  of  a  computer  program  for  tire 
performance  calculations  based  on  the 
mathematical  model 

0  Verification  of  the  model  bv  comparisons 
of  predicted  tire  performance  with  avail¬ 
able  experimental  data. 
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IDEAL  FORMULATION  OF  TIRE -SOIL  INTERACTION 


I  I . 


Theor  e  i  i  ca  1 1  v  rigorous  Simula  Lion  of  physical  phenomena  is 
i  llusorv  since  the  characterization  oi  the  physical  properties  of 
materials  and  the  description  of  the  natural  law.  that  govern 
their  behavior  are  themselves  idealizations.  Rigorous  solutions 
exist  onLy  for  hypothetic  t 1  materials  and  conditions.  Thus,  in 
principle,  the  physical  phenomenon  of  tire-soil  interaction  can¬ 
not  tie  treated  in  a  theoretically  rigorous  manner.  An  ideal 
formulation  would  be  one  that  would  yield  a  rigorous  solution 
for  the  hypothetical  case  that  both  tire  and  soil  behave  as  Lheir 
material,  characterization  postulates.  For  such  a  hypothetical 
condition  it  would  x  possible  to  write  a  symbolic  formulation 
for  the  problem  on  tne  basis  of  the  equations  of  motions  and  the 
boundary  conditions  at  the  interface.  A  general  solution  of  this 
problem,  however,  would  present  a  formidable  task  and  may  he  re¬ 
garded  only  as  the  ultimate  goal.  Solutions  of  special  cases 
that  are  considered  attainable  with  present  state  of  the  art 
methods  involve  the  following  assumptions: 

The  tve  travels  in  a  straight  line  at  a 
constant,  low  velocity 

0  A  "steady  state"  exists  in  the  soil 

*  The  tire  load  is  constant  —  that  is,  there 

is  no  interaction  between  tire  and  suspension 
sys tern 

c  The  terrain  is  even, 

To  formulate  the  tire-soil  interaction  for  these  conditions, 
a  tire  model  is  needed  for  the  ca leu  la ti on  of  tire  deformation 
under  a  system  of  applied  stresses.  Finite  element  models  of  tires 
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that  have  been  developed  for  other  purposes  would  generally  re¬ 
quire  extensive  further  development  if  they  were  to  be  applied  to 
die  problem  of  Lire-.->oil  interaction.  While  it  •  ould  be  desirable 
to  have  such  a  model  for  the  study  of  tire  deformation  under  vari¬ 
ous  conditions  occurring  in  tire-soil  interaction  situations, 
finite  element  models  generally  require  large  computer  capacity 
.and  proportionate  Lime  for  the  computation  of  deformations  under 
a  single  set  of  conditions.  Since  in  a  predicLive  tire-soil 
model  such  computations  would  he  needed  repeatedly  because  of 
the  iterative  nature  of  the  prediction  technique,  the  use  of  finite 
element  tire  models  was  deemed  undesirable.  Theoretical  considera¬ 
tions  as  well  as  studies  cf  tire-soil  interaction  experiments  indi¬ 
cated  that  for  the  purpose  of  performance  prediction  a  refined  finite 
element  tire  model  was  not  necessary,  and  that  a  simple  tire  model 
responsive  to  changes  in  soil  stiffness  could  be  successfully  used. 
The  experimental  information  available  on  the  geometry  of  tires  in 
various  soil  conditions  and  the  development  of  the  approximate  tire 
model  arc  described  in  subsequent  sections. 
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Ill . 


KX PER I MENTAL  INFORMATION  ON  TIRE -SOIL  INTERACTION 


Of  the  numerous  tiro  tests  that  have  been  performed  by  various 
researchers  over  the  years,  the  most  valuable  for  the  development  of 
tire  model  are  those  where  interface  stresses  and  Lire  deforma¬ 
tions  oi  tires  moving  in  soils  were  directly  measured.  The  experi¬ 
mental  data  on  the  deflection  of  tires  moving  in  soft  soils  and  the 
distribuiion  of  interface  stresses  obtained  in  various  WES  research 
project^  (Refs.  1  through  6)  are  the  most  instructive  available  in¬ 
formation  on  the  role  of  Lire  deflection  in  tire-soil  interaction 
and  served  as  a  basis  for  the  development  of  the  tire-soil  model 
presented  in  Section  VII. 

Interface  stress  and  tire  deflection  measurements  performed  by 
others  and  available  in  the  Literature  (Refs.  7  'Trough  10)  were 
also  studied  and  used  in  the  development  of  the  tire-soil  model. 

Direct  measurements  of  Lire  deflection  and  interface  stresses 
are  expensive  anti  time  consuming,  and,  therefore,  only  few  were 
pe  r  l  tinned  .  On  the  other  hand,  performance  tests,  without  direct 
geometry  and  stress  measurements,  have  be  n  performed  for  a  wide 
variety  of  conditions.  Since  Lire  geometry  changes  wiLh  the  slip 
rate,  those  tests  where  performance  was  measured  over  a  range  of 
^lip  offered  tine  experimental  information  necessary  for  the  develop¬ 
ment  and  validation  of  the  tire-soil  model.  At  our  request,  WES 
furnished  data  on  such  tests  (Ref.  11)  on  various  types  of  tires. 

A  summary  of  tests  performed  at  WES  in  Yuma  sand  is  given  in  Table  1 
anti  of  those  performed  in  Buckshot  clay  in  Table  2.  Results  of 
these  tests  are  shown  graphically  in  conjunction  with  the  evalua¬ 
tion  of  predictions  by  the  tire-soil  model.  More  detailed  informa¬ 
tion  on  the  performance  of  these  tests  is  contained  in  WES  Technical 
ReporL  No.  3-fS6b  (Ref.  12). 
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The  WIN  tes  t  s  were  performed  in  either  Liu  purely  frioLional 
Yuma  and  or  in  the  purely  cohesive  Buckshot  lay.  Ln  order  to 
hev>  some  information  on  tire  performance  in  f ric  t  i  ona  i-cohesi.  ve 
soi 1  ,  results  of  tire  tests  performed  in  two  clayey  soils  at  the 

National  Tillage  Machinery  Laboratory,  Auburn,  Alabama,  were  ob¬ 
tained  (Table  !)•  These  tests  were  performed  for  comparative 
studies  rathe-''  than  performance  prediction,  and,  therefore,  evalua¬ 
tion  ot  oit  strength  had  to  be  done  in  retrospect  from  Lhe  avail¬ 
able  data  on  soil  proper  ties . 
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HE NEKA E  APPROACH  TO  THE  DEVELOPMENT  01’  TIRE -SOLE  MODE  1 


r  v . 


In  view  of  the  paucity  of  specific  deflection  and  stress 
measurements  and  the  wide  variety  of  tire  performance  tests, 
the  following  approach  tc  the  development  of  tire-soil  model  was 
pursued . 

<•  Establish  qualitative  relationships  between  Lire 
deflections  and  soil  properties  and  between  in¬ 
terface  stresses  and  inflation  pressure 

'*  Modify  the  validated  rigid  wheel -soil  model 

(Ref.  13)  so  as  to  allow  for  the  influence  of 
lire  deflection  and  inflation  pressure.  In  the 
mode  1  allow  for  p  sufficient  number  of  free 
parameter,  so  that  existing  deflection  and 
stress  measurements  may  be  matched 

J  Analyze  the  effect  of  free  parameters  and  other 
assumptions  on  tire  performance,  modify  assump- 
1  ions,  and  establish  values  for  free  parameters 
on  the  basis  of  comparisons  with  tire  performance 

test  s 

Establish  methods  to  determine  model  parameters 
for  prediction  purposes. 
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V  . 


TIRE  DEFLECTION  IN  SOFT  SOILS 


The  tire  deflection  measurements  carried  out  at  WES  under 
various  loading  and  soil  conditions  and  reporteu  in  Refs.  1 
through  6  have  conclusively  shown  that  deflection  of  the  cross 
section  of  tires  as  well  as  that  of  their  centerline  in  the  plane 
of  motion  is  dependent  on  both  the  stiffness  of  the  tire  and  the 
stiffness  of  the  soil.  The  deflection  of  the  cross  section  of 
the  tire  influences  tire  performance  primarily  by  enlarging  the 
width  of  the  contact  area.  This  width  change  is  within  a  narrow 
range,  its  lower  limit  being  the  width  of  the  unloaded  tire  and 
its  upper  limit  the  width  of  the  tire  loaded  on  a  rigid  surface. 
These  limiting  widths  are  generally  available  or  can  be  easily 
measured;  the  width  of  a  tire  loaded  in  soft  soil  may  be  estimated 
from  these  limits  with  sufficient  accuracy.  Since  the  proposed 
tire-soil  model  is  a  two  dimensional  one,  effects  of  transverse 
deflection  are  considered  only  by  changes  in  width. 

The  deflection  of  the  centerline  of  the  tire  in  the  plane  ot 
motion  determines  the  average  inclination  of  the  contact  area  to 
the  horizontal  and  thereby  influences  tire  performance  profoundly. 

The  inclination  of  the  contact  area  acts  as  if  the  tire  were 
to  climb  a  slope.  For  a  rigid  wheel  the  average  inclination  of 
the  contact  area  is  the  angle  bisecting  the  entry  and  exit  angles 
and  may  easily  be  in  the  range  of  10  to  20  degrees.  On  the 
other  hand ,  a  very  flexible  tire  may  flatten  out  appreciably  so 
that  the  contact  area  becomes  almost  horizontal.  The  inclination 
of  the  contact  area  affects  the  pull  coefficient  directly,  approxi 
mately  by  the  value  of  the  tangent  of  the  angle  of  the  inclination 
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The  above  considerations  show  the  imporlance  of  cenLei line 
'(.oiiKlry  in  a  lire-soii  model.  The  experiments  performed  at  WES 
show  the  variation  of  centerline  geometry  with  slip,  soil  strength, 
ind  inflation  pressure.  From  these  experiments  the  following 
qualitative  conclusions  were  drawn. 

The  centerline  geometry  of  a  tire  is  affecLed  by  slip,  infla¬ 
tion  pressure,  and  soil  strength.  The  combined  effect  of  tire  and 
soil  stiffness  may  be  graphically  represented  as  shown  in  Fig.  1. 
Tire  deformation  is  negligible,  and  the  tire  behaves  as  a  rigid 
wheel  if  the  tire  stiffness  is  great  relative  to  Lhat  of  the  ground. 
In  the  other  extreme,  tire  deformation  is  comparable  Lo  Lhat  on  a 
rigid  surface  if  the  stiffness  of  the  soil  reiacive  Lo  that  of  the 
tire  is  great. 

Another  finding  related  to  the  geometry  of  centerline  is  that 
the  resultant  of  normal  stresses  al  the  Lire  soil  interface  nasses 
through  the  centerline  oi  the  whee 1  axle  (Ref.  A). 

These  qualitative  relations  defining  the  effect  of  tire  and 
soil  stiffness  on  centerline  geometry  were  the  basis  for  the  de¬ 
velopment  of  the  tire-soil  model. 
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VI  . 


STRESSES  BENEATH  A  '1  IRE  MOVING  IN  SOFT  SOIL 


Tire  deflection  aflecLs  tire-soil  interaction  not  only  by  its 
i  ilcct  on  the-  geometry  ol  Lite  contact  area,  hut  also  by  relieving 
the  stresses  that  would  develop  in  the  soi  1  if  the  Interlace  were1 
unde  i  orr.iab  1  e  ,  Stress  measurements  at  the  tire-soil  interlace  con¬ 
firm  this  latter  effect  and  give  an  indication  of  the  magnitude  of 
s tress  relief. 

Tire-soil  interface  stress  Measurements  at  WES  were  first 
made  on  unyielding  surfaces  with  sensor  placement  in  the  unyield¬ 
ing  surface  (Ref.  5).  Results  of  these  measurements  are  of  in¬ 
terest  for  Lire-soil  interaction  studies  because  stresses  measured 
on  an  Unyielding,  surface  represent  the  upper  limit  of  stresses 
that  would  develop  in  a  soil  that  yields  relatively  little  under 
lire  load.  The  general  pattern  of  stress  distribution  observed 
in  'hose  tests  showed  a  center  portion  in  Lhe  contact  area  with 
fairly  uniform  stress  distribution  and  stress  oncen  tra  t  ions  , 
called  "edge  stresses,"  at  the  perimeter  of  the  contact  area. 

These  edge  stresses  are  related  to  the  sidewall  stiffness  of  the 
tire,  while  the  magnitude  of  the  average  center  portion  stresses 
is  related  to  Lhe  inflation  pressure  of  the  tire. 

Measurements  of  interface  stresses  in  soils  were  also  made 
by  VandenBerg  and  Gill  (Ref.  h) .  These  measurements,  using  smooth 
tires,  indicate  a  stress  distribution  pattern  similar  to  that  ob¬ 
served  on  unyielding  surfaces.  The  magnitude  of  the  uniform  pres¬ 
sure  in  the  center  portion  depends  on  the  inflation  pressure  and 
is  generally  somewhat  higher. 

Freitag  et  al.  (Ref.  4)  investigated  the  distribution  of  nor¬ 
mal  stresses  in  the  contact  area  of  both  towed  and  powered  tires 
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mikiud  to  various  pressures.  Tests,  carried  out 
' 11111  l~  * 0 v  '  i  nt*  i (,<J Led  tint  the  maximum  normal  stress 
exceeded  only  sjj  ht  ly  the  inflation  pressure. 
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VII.  CONCEPT  OF  TIRE -SOIL  MODEL 


The  de f oriua *■  ion  and  stress  measurements  discussed  briefly  in 
Section  VI  indicate  the  complexity  of  the  tire-soil  interaction 
problem.  The  shape  of  the  Lire  and  geometry  of  the  contact  area 
depend  not  only  on  the  properties  of  the  tire,  but  also  on  the 
properties  of  soil  and  on  the  loads  applied.  The  stresses  mea¬ 
sured  in  the  contact  area  are  far  from  uniform;  stress  concentra¬ 
tions  occur  at  the  edges  of  the  contact  area.  Obviously,  all 
these  variations  cannot  be  considered  in  any  workable  t.ire-soil 
model,  and  simplifications  are  required.  An  appropriately  sim¬ 
plified  model  often  yields  sufficiently  accurate  results,  as  many 
computational,  methods  in  engineering  demonstrate. 

To  decide  what  simplifications  can  be  undertaken  in  a  model 
without  jeopardizing  its  accuracy  and  usefulness,  it  is  expedient 
to  consider  the  Lire  as  a  free  body  and  to  assess  the  effect  of 
possible  simplifications  on  the  performance  of  the  Lire.  The  edge 
stresses  in  the  contact  area,  as  experiments  indicate,  are  sym¬ 
metrical  both  crosswise  and  lengthwise.  This  symmetry  allows  one 
to  consider  average  stresses  across  the  tire  width  without  any  sig¬ 
nificant  loss  of  accuracy.  Likewise,  edge  stresses  may  be  smoothed 
lengthwise  and  the  resulting  torque,  load,  and  drawbar  pull  still 
may  be  reasonably  close  to  the  actual  values.  Conversely,  it  is 
important,  to  duplicate  the  defLected  shape  of  the  tire  and  its 
orientation  to  the  ground  surface  as  closely  as  possible.  In  the 
summation  of  the  interface  stresses  for  the  computation  of  drawbar 
pull,  the  inclination  of  int'rface  elements  relative  to  the  ground 
surface  cannot  be  neglected.  Depending  on  the  inclination  of  the 
element,  normal  stresses  on  the  interface  yield  a  component  (plus 
or  minus)  in  tht  direction  of  the  drawbar  pull  that  may  or  may  not 
be  significant  relative  to  the  component  of  shear  stresses. 


Tin-  Lin.  -  :  l  modi  i,  onecived  on  Lhe  bam  ol  these  con- 

s  i  dura  t  i  ons  ,  is  shown  in  !•  i .  2.  The  ossonli.il  loatuios  ol  i  ho 

;  ,o(lu  i  an-  as  I  o  I  1  ova  :  The  t  »  i*i-  is  assumed  Lo  have  a  ooiisl.ini  width 
in  bolh  Lhu  undo  1  onurd  and  do  formed  s  La  Los.  The  stresses  amiss 
Lhe  Lire  width  are  assumed  Lo  lx-  uniform  so  that  Lhe  lirc-soiJ  in- 
ieruction  problem  can  he  treated  as  Leo  dimensional.  Tire  deforma¬ 
tion  is  represented  by  Lhe  shape  of  Lhe  Lire  in  the  center  plane  in 
lhe  direction  of  travel.  Tire  shape  is  assumed  Lo  be  lhu  same  in 
all  parallel  planes  The  centerline  geometry  is  assumed  to  cons  is L 
of  lv.o  curvilinear  segments  separated  by  a  linear  or  flat  section. 

[L  is  assumed  that  the  Lire  sLarts  to  deform  ahead  of  the  entry 
ingle  and  reaches  its  original  form  pasl  Lhe  exit  angle  ■  , 

:  ^  shown  in  Fig,.  2.  The  curvilineai  segments  are  logarithmic 
pi  ra  1  s  vjlth  Lhe  radii  decreasing  according  Lo  the  following  re¬ 
in!  ions  hip 

r  =  R  e  (  °}  (1) 

v.  here 

li  =  radius  of  unloaded  Lire 
r  =  deflected  radius 
=  constant 
=  central  angle 

=  initial  central  angle 

o 

The  effect  of  tire  deflection  on  soil  reaction  is  represented 
in  the  proposed  tire-soil  model  in  two  ways:  Soil  failure  condi¬ 
tions  as  determined  by  plasticity  theory  govern  the  interface  nor¬ 
mal  stresses  in  a  front  and  rear  failure  zone.  Deflection  of  the 
centerline  affects  the  magnitude  of  these  stresses;  this  effect  is 
taken  into  account  by  using  the  deflected  geometry  as  boundary 
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ondilion  in  tin  numori  al  computation.  FurLhe nnore  ,  in  the  model 
is  assumed  LhuL  Lire  deflection  limits  L he  normal  stresses  that 
au v  arise  a !  In  inte’  'ace.  This  idea  was  originally  proposed  by 
kkei  hi  in  i  for  a  towed  t.  ire  model  (Ref.  Id).  The  may, ni  Lucie 

■  : i :  1;  .a  i  iny,  pressiuc-  depends  on  the  inflation  pressure  and 

ar  s'  i::ness.  In  the  proposed  Li  re-soil  model  the-  limiting 

,  ressets  i  nnsidt  red  .is  a  ree  parameter  whose  value  is  to  be 
t  si  ubl  i  shed  by  . .  pa  rattle  tri  n  analysis  of  exper  i  morn  a  1  results.  The 
toll  oi  hear  st  tosses  generated  by  the  applied  torque  is  implicit 
in  the  del  in  that  the  ge  rv  ol  the  slip  line  fields  and  the 
ussivi  iteci  nterLaee  norma  resses  depend  on  them.  The  normal 
stress  of  the  inner  end  the  rear  failure  zone  equals  the  limiL- 
i  n .  pres  ire  p  j  .  If  t  he  rise  oi  the  normal  stresses  in  this  zone 
is  mode:  <  r  ;  1  ov-  be  muse  of  the  app  li  ed  shear  stresses ,  the 

:  i  iure  ne  ■  :  end:  .  urther  t  award  the  front  zone,  and  the  zone 
;  uni  for:  norma  !  s  t  re.sses  where  the  soil  is  not  in  failure  con¬ 
dition  is  rodueeu  or  entirely  eliminated.  In  such  cases  the  tire- 
soil  .ode  1  Ik  turn  tilar  to  the  rigid  wheel-soil  i  oelel  . 

lht  key  mature  ol  the  proposed  tire-soil  model  is  the  deter¬ 
mination  of  the  eons'  mL  in  Eq.  1).  This  constant  determines 

the  curvature  oi  the  I  ogariLhmi  c  spiral  and  is  a  measure  of  tire 
defie  lion.  As  discussed  in  Section  V,  tire  deflection  in  soft 
soils  depends  on  the  relative  stiffness  of  the  tire  to  that  of  the 
soil.  Tire  stiffness  depends  on  the  inflation  pressure,  while  soil 
stiffness  is  directly  related  to  the  strength  of  the  soil.  The 
coefficient  in  Eq  (1)  should  reflect  the  combination  effect 

of  tire  and  soil  stiffness.  In  the  mathematical  formulation  of 
the  model  this  is  accomplished  in  the  following  way ;  The  arc  length 
of  the  interface  of  the  forward  failure  zone  depends  on  bolh  the 
limiting  pressure  and  the  soil  strength.  In  strong  soils  the  rise 
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of  the  normal  stresses  is  steep,  and  the  limiting  pressure  is 
reached  over  a  short  arc  length;  whereas  in  a  weak  soil  a  longer 
arc  lengLh  is  needed  to  reach  the  same  limiting,  pressure.  II 
soil  conditions  are  the  same,  the  arc  lengLh  depends  on  l  lie  limit¬ 
ing  pressure;  low  limiting  pressure  will  result  in  small  are  length, 
while  high  limiting  pressure  will  result  in  long  arc  length.  Thus 
the  arc  length  of  ihe  forward  failure  zone  is  governed,  at  least 
qualitatively,  by  the  same  interrelationship  lhat  governs  Lire  de¬ 
flection  according  to  the  experiments  illustrated  schematically  in 
Fig.  1.  In  order  to  introduce  this  interrelationship  into  the 
mathematical  formulation  of  the  proposed  tire-soil  model,  a  de¬ 
flection  coefficient,  c,  is  introduced  that  defines  the  shorten¬ 
ing  of  the  tire  radius  by  deflection  at  a  specified  central  angle, 

(a  ) ,  as  fol lows : 
c 

r  -  eR  .  (2) 

The  coefficient  is  calculated  for  any  given  c  value  from 

Eq.  (1).  By  this  method  of  calculation,  depends  on  the  arc 

length  of  the  front  failure  zone  and  reflects,  at  least  qualitatively 
the  desired  interrelationship  between  Lire  deflection  and  soil  stiff¬ 
ness.  In  the  development  stage  the  angle  •  Lo  which  Eq .  (2)  re- 

l 

fers  was  considered  as  a  free  parameter.  The  analysis  of  Lire  per¬ 
formance  experiments  showed  that  good  simulation  can  be  achieved  if 
=0.3  ,  is  assumed . 

C-  (i 
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VIII.  DEVELOPMENT  OF  MATHEMATICAL  TIRE-SOIL  MODEL 


The  concept  of  the  tire-soil  model  described  in  Section  VII 
takes  into  account,  at  least  qualitatively,  the  interaction  ef¬ 
fects  between  Lire  and  soil  evidenced  in  various  experiments.  In 
the  mathematical  model  it  was  necessary  to  express  the  interrela¬ 
tionships  quantitatively  and  check  the  validity  of  these  relation¬ 
ships  against  experimental  data.  The  first  step  in  this  direction 
was  to  write  a  computer  program  for  the  model  and  check  whether  the 
model  was  capable  of  duplicating  experimentally  determined  center- 
line  geometries  and  stress  distributions.  In  this  matching  pro¬ 
cedure,  the  as  yet  undefined  free  parameters  were  selected  by  trial 
and  error.  While  the  assessment  of  whether  an  acceptable  duplica¬ 
tion  of  the  experimental  data  was  achieved  was  somewhat  subjective, 
this  trial  and  error  exercise  was  also  useful  in  obtaining  a  feel 
for  the  effect  of  the  variables  on  the  tire  geometry  and  stresses. 
In  many  cases  it  was  fo und  Lhat  close  duplication  could  he  achieved 
by  various  combinations  of  the  free  parameters,  pointing  to  the 
possibility  of  reducing  the  number  of  free  parameters  in  the  mathe¬ 
matical  model.  The  basic  elements  in  the  computer  program  used  in 
this  development,  such  as  the  subroutine  for  the  calculation  of 
slip  line  ields,  were  essentially  the  same  as  in  the  final  program 
for  the  pc  rfornunce  prediction  presented  in  Appendix  B,  but  many  of 
the  iteration  procedures  were  operated  with  manual  interaction. 
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Although  the  computer  program  in  the  above  development  stage 
with  its  capability  to  duplicate  experimental  data  was  valuable  in 
itself  for  the  analysis  of  tire-soil  interaction,  for  the  purpose 
of  tire  performance  prediction  it  was  necessary  to  establish  its 
general  validity  and  develop  relationships  among  the  free  and  input 
parameters  and  methods  to  determine  the  input  parameters  from  available 
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information  on  tiro  and  soil  for  t ho  prod io t  i on  ol  po r f ormanoo . 
Since  experimental  information  on  tiro-soil  inloraction  that 
covered  a  u  i  d  o  ran^o  of  conditions  and  a  number  of  tiro  typos  uas 
liniLod  to  tiro  porlormanco  data,  a  method  was  dovisod  to  utili/.o 
•_ h o s o  data  for  these  purposes.  The  computer  program  at  this 
stage  established  a  multivariate  relationship  between  tiro  per¬ 
formance  and  a  number  of  free  and  input  parameters.  The  problem 
uas  to  find  values  for  those  parameters  so  that  »  he  experimentally 
determined  performance  (.lata  be  matched  for  all  conditions  tested. 
The  available  experimental  data  lisled  in  Tables  1  and  2  repre¬ 
sented  3.'  different  conditions;  each  of  the  38  tests  '’onlained 
load,  pull,  torque,  and  sinkage  data  at  several  values  of  slip.  The 
•orrpuLer  Lime  required  to  compute  the  performance  parameters  for 
the-  possible  combinations  of  free  and  input,  parameters  for  each 
slip  rate  in  the  3  tests  on  the  PDF- 10  computer  was  economically 
prohibitive,  and  recourse  was  made  to  the  use  of  minicomputers  for 
this  purpose.  In  an  independent  research  project  of  the  Research 
Department  methods  were  developed  to  run  relatively  large  programs 
on  minicomputers  b /  storing  segmenLs  of  the  program  on  disks  and 
bring,  selected  portions  into  the  core  when  needed  (Ref.  15).  As  a 
result  of  this  operation  and  the  lower  computation  speed  available 
with  minicomputers  the  running  time  is  about  two  to  six  times 
longer  than  on  the  PDP-10  or  similar  computers,  but  the  cost  of 
computation  is  minimal.  The  program  was  easily  adapted  to  the  re¬ 
search  department's  NOVA  800  minicomputer;  the  main  program  and  the 
subroutine  constituted  two  segments  that  the  core  was  capable  of 
accommodating.  The  longer  running  time  was  of  no  disadvantage  in¬ 
asmuch  as  loops  for  the  desired  variations  of  parameters  were  in¬ 
corporated  into  the  program,  and  the  extended  program  was  run  over¬ 
night  without  any  supervision.  This  way  performance  calculations 
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wert’  made  lor  over  IS, 000  combinations  of  the  parairn’tors  without 
incurring  costs  oLher  than  minimal  depreciation  and  maintenance. 

Even  though  the  use  of  the  minicomputer  r  ide  an  extensive 
analysis  of  the  program  feasible,  it  was  desirable  Lo  reduce  the 
number  of  combinations  of  the  various  parameters  as  much  as  possi¬ 
ble.  To  this  end  it  vas  decided  to  a  1  low  variation  of  certain 
parameters  only  if  no  combination  of  the  ot  ler  parameters  yielded 
satisfactory  results.  This  restriction  of  the  variation  of  the 
parameters  involved  the  following  assumptions: 

The  •'  and  angles  (Fig.  2)  were  assumed 

to  be  equal 

The  curvature  of  the  logarithmic  spiral  section 
was  assumed  to  be  the  same  in  the  front  and  rear 
sec t ion 

-  It  was  assumed  that  the  soil  strength  parameters 
■;  and  e  were  uniquely  determined  by  the  cone 
index  or  its  gradient  in  the  experiments  lisLed 
in  Tables  1  and  2.  The  following  relationships 
between  strength  and  cone  index  parameters  were 
es  tab fished : 


Yuma  sand:  J  (%)  =  71.1  log  (CGR"  +  11.33 

(3) 

T  =  arc  tan^l/(1.64  -  0.6H  D  ) 


Buckshot  c lay : 

C(psi)  =  Cl/ 12. 5 
<1  (°  )  =  Cl/  4 


(4) 


.  v.ti?  *v.O i  ■  r- •  £&<!*‘eVb. » 
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Since  Lhe  es  Lab  i  is  hu.en  L  of  lIicsc  relationships  ls  closely 
related  to  the  genera  1  problem  of  the  determination  of  strength 
properties  ol  soils  for  mobility  purposes,  iL  is  discussed  in 
detail  under  th.it  heading  in  Appendix  A. 

Tiie  effect  of  model  parameters  on  performance  simulation  was 
analyzed  by  a  systematic  variation  of  the  following  parameters: 

0  Deflection  coefficient  <■ 

°  Limiting  pressure  p. 

■-  Any  le  "  = 

-  Angle 

t 

'■  Initial  values  of  a. no 

r  d 

In  these  performance  simulation  analyses  the  main  objective 
v.as  to  simulate  the  puLl  coc  f f icient -slip  relationship.  Simula- 
t  ion  of  sinkage  ..as  considered  only  as  a  secondary  objective, 
partly  because  sinkage  measurements  arc-  somewhat  an  ertain  (sink- 
age  was  not  directly  measured  in  the  tests  but  derived  by  an  ex¬ 
perimental  formula  from  the  vertical  movement  of  the  hub).  Com¬ 
parison  of  computed  and  measured  sinkage  was,  however,  very  useful 
In  deciding  wh'lher  discrepancies  in  pull  coefficients  occurred 
because  of  had  centerline  geometry  simulation  or  for  other  reasons. 
Generally  speaking,  good  simulation  of  the  pull  coefficient  was  not 
possible  without  fairly  good  sinkage  approximation.  Computed  and 
measured  torque  coefficients  generally  agreed  well  when  the  pull 
coefficient  simulation  was  good;  reasons  for  exceptions  are  dis¬ 
cussed  in  Section  X. 

The  performance  simulation  analyses  resulted  in  the  establish¬ 
ment  of  initial  values  of  =  1CT  and  ■  .  =  f  . 

r  d 
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Ln  siron”  cohesive  soils,  I  loses  its  meaning  when  the 
strength  of  the  soil  is  such  that  there  is  no  failure  condition 
in  Lhe  front  field  anci  only  a  rear  plastic  failure  zone  develops 
the  angle  of  separation  of  the  front  field.  I  i  such  a  case,  good 
simulation  results  if  the  entry  angle  is  assumed  as  one  and  a 
half  times  Lhe  rear  angle. 

For  the  value  of  the  angle  •  the  relationship  >  =  '’t|/2 

was  found  to  result  in  good  simulation. 

The  values  established  above  for  these  angles  are  not  neces¬ 
sarily  the  best  ones,  but  were  found  to  yield  acceptable  perfor¬ 
mance  simulation  in  the  cases  analyzed. 

In  the  concept  of  the  tire-soil  model,  Lhe  angle  of  the  de¬ 
veloped  interface  friction,  £,  is  the  independent  variable  that 
enters  into  the  computation  of  slip  line  field.  The  angle  5  is 
related  to  slip  by  the  following: 

-(j+j0)/K 

tan  d  =  tan  5  i  1  -  e  (5) 

maxV 

This  equation  is  the  same  type  as  proposed  by  Janosi  and 
Hanamoto  for  the  relationship  between  mobilized  shear  stress  and 
slip  for  tracked  vehicles  (Ref,  16).  The  slip-shear  parameters 
iii  iq,  ( 5 )  we  re  not  known  beforehand  and  had  to  be  assumed  so  that 
a  slip  value  corresponding  to  the  f  value  could  be  determined 
and  the  computed  pull  coefficients  compared  with  the  measured  ones. 

With  the  values  of  the  various  angles  established  by  Lhe 

analysis  for  model  development,  a  systematic  analysis  of  each  of 
the  31  tire  performance  tests  was  made  to  determine  Lhe  values  of 
the  parameters  e,  p^  ,  j  ,  and  K  that  yielded  good  simulation. 
Figures  3  through  8  are  examples  of  good  simulation  of  pull  coeffi¬ 
cients  obtained  In  this  systematic  analysis.  The  excellent  agreement 


29 


0.8 


Pull 

Coef  f 


’0.2 


~r  Computed 
+  Variable  Slip  Test 


Tire:  9.00-14  Load:  450  lbs  Infl.  Pressures  16.4  psi 

Defl:  157.  e  -  0.97  P/  -  16.5  psi 

Cone  Index  Gradient:  8.9  Friction  Angle:  42.4° 

Slip  Parameters:  i  ■=  0.47  K  =  0.184 

'  o 

Fig.  3  Simulation  of  Pull  Performance  with  Freely 
Selected  Parameters,  Test  No.  1 
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*  Computed 
4-  Variable  Slip  Test 


o  o 

-  ♦ 

♦  ♦ 

Pull  ♦ 

Coeff  ♦ 


Tire:  6.00-16  Load:  455  lbs  Infl.  Pressure;  10.3  psi 

Def  1 :  257.  e  =  0.88  ??  ■=  10  psi 

Cone  Index  Gradient:  15.6  Friction  Angle:  45.4° 

Slip  Parameters:  j  =0.03  K  =  0.07 

r  j  q 

Fig.  4  Simulation  of  Pull  Performance  with  Frcely 
Selected  Parameters,  Test  No.  4 
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*  Computed 
+  Variable  Slip  Test 


I 


[■ 


Pull 

Coe  f  f 
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0  Slip,  7,  80 

Tire:  9.00-14  Load:  890  lbs  Infl.  Pressure:  16.4  osi 

Defl:  2  57,  e  =  0.88  P  ^  ■=  16.0  psi 

Cone  Index  Gradient:  13.8  Friction  Angle:  44.7° 

Slip  Parameters:  j  =  0.03  K  =  0.12 
r  J  o 

Fig.  5  Simulation  of  Pull  Performance  with  Freely 
Selected  Parameters,  Test  No.  9 
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w  Computed 
+  Variable  Slip  Test 


Pull 

Coeff 


Slip,  %  80 

oea;  P/30.1SspslInfl-  Pressure:  7'50pSi 

sup  sss-j  -xo-'m04  Trsv*1"  — 


Fig.  6  Simulation  of  Pull  Performance  with  Freely 
Selectee  Parameters,  Test  No.  11 
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Slip,  7o 


I 

80 


Tire:  4.00-20  Load:  340 
Defl:  3  5%  e  =  0.92  p.  - 
Cone  Index  Gradient:  25.6 
Slil  Parameters:  j  =  0.03 


lbs  Infl.  Pressure:  11  psi 
10  psi 

Friction  Angle:  46  2° 

K  =  0.18 


Fig.  7  Simulation  of  Pull  Performance  with  Freely 
Selected  Parameters,  Test  No.  22 
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•'  Computed 
+  Variable  Slip  Test 
x  Constant  Slip  Test 


Pull 
Coef  E 


♦ 


4 


0 
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0 


4 


-0.2 


-20 


Slip,  7o 


SO 


Tire:  9.00-14  Load:  455  lbs  Infl.  Pressure:  7.5  psi 

Defl :  2  57,  e  =  0.87  P;  =  8.5  psi 

Cone  Index  Gradient:  11.8  Friction  Angle  =  43.8° 

Slip  Parameters:  jG  =  0.04  K  =  0.11 


Fig.  8  Simulation  of  Pull  Performance  with  Freely 
Selected  Parameters,  Test  No.  25 
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between  experimental  and  calculated  values  was  indicative  of  the 
capability  of  the  model  to  simulate  tire-soil  interaction  and  in 
most  cases  was  within  the  probable  limits  of  accuracy  of  Lhe  ex¬ 
periments.  Obviously,  there  was  no  need  for  further  refinement 
of  the  model,  hut  relationships  between  parameter  e,  p^ ,  j  , 
and  K  and  measurable  Lire  and  soil  properties  still  had  to  be 
developed  for  prediction  purposes.  An  extensive  study  of  Lhe 
freely  selected  parameters  that  yielded  good  simulation  was  made 
and  Lhe  fo 1  lowing  rules  Lor  the  estimation  of  these  parameters 
were  developed. 

Estimation  of  the  Deflection  Coefficient  c 

The  deflection  coefficient  c  is  closely  related  to  the  de 
sign  deflection  coefficient  2f/d  used  by  WES  and  can  be  deter¬ 
mined  by  a  deflection  measurement  on  a  rigid  surface.  Table  4 
contains  the  relationships  for  the  estimation  of  coefficient  e 
for  the  various  sizes  of  tires  in  terms  of  v:  =  1  -  2 6 / d . 

TABLE  4.  ESTIMATION  OF  THE  DEFLECTION  COEFFICIENT’ 

«  FOR  VARIOUS  SIZE  TIRES 


Tire  Size 

Relationship  for  Estimation,  e  = 

o.00-14 

1.24x  -  0.195 

6.00-16 

1.62k  -  0.59 

4.00-7 

1  .  1 4 >:  -  0.115 

4.00-20.0 

0 . 91k  +  0.085 

31  x  15-13 

0.98k  +  0.065 

The  relationships  listed  in  Table  4  for  the  estimation  of  a 
are  shown  in  Fig.  9.  It  can  he  seen  that  they  are  very  similar  to 
each  other  and  allow  an  approximate  estimate  of  <r  for  tire  sizes 
other  than  those  listed  in  Table  4.  If  no  other  information  is 
available,  the  relationship  r  =  k  shown  by  a  heavy  dashed  line 
in  Fig.  9  can  serve  as  a  first  approximation. 

Estimation  of  the  Limiting  Pressure,  p^ 

It  was  found  that  the  limiting  pressure  p^  may  be  estimated 
on  the  basis  of  inflation  pressure  for  all  tire  sizes  by  the  fol¬ 
lowing  relationship: 

p^(psi)  =  0.64  x  p  +  4  (6) 

The  Pj  values  determined  by  the  above  relationship  are 
higher  than  p^  up  to  about  11  psi  and  lower  than  p.  above 
11  psi.  (The  Pj  values  should  noL  be  interpreted  as  the  average 
pressure  in  a  rigid  surface.)  The  p^  value  is  a  model  parameter 
that  refers  to  the  limiting  pressure  that  yields  acceptable  tire- 
soil  interaction  simulation  with  the  proposed  tire-soil  model.  It 
may  be  interesting  to  note,  however ,  that  some  tires  at  high  values 
of  inflation  pressure  exhibit  average  pressures  lower  than  the  in¬ 
flation  pressure  (Ref.  17). 

Estimation  of  the  Slip-Shear  Parameters  i  and  K 

_ _ _ _ _o _ 

The  equation  for  the  relationship  between  mobilized  shear 
strength  and  slip  proposed  by  Janosi  and  Hanamoto  (Ref.  16)  was 
based  on  the  analogy  of  tracks  and  direct  shear  tests  from  which 
the  parameter  K  in  the  proposed  equation  could  be  evaluated.  In 
applying  the  same  concept  to  wheels  it  was  necessary  to  include  the 
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inici  value  ]  ^  in  the  slip  term  since  experiments  showed  a 
negative  slip  at  the  towed  point  where  the  developed  shear  stress 
is  zero.  In  the  case  of  wheels  and  tires  the  analogy  with  the 
direct  shear  test  is  not  evident  since  it  is  difficult  to  define 
aid  measure  the  differential  dis  placement  between  points  on  the 
wheel  surface  and  soil.  Theoretical  modeling  of  the  slip-shear 
phenomenon  in  both  shear  strength  tests  and  tire-soil  interaction 
would  be  necessary  to  develop  a  theoreLica1  basis  for  tK  estima¬ 
tion  of  throe  parameters.  In  view  of  the  lack  ol  such  theory  ex¬ 
tensive  studies  were  made  to  determine  whether  there  is  a  correla¬ 
tion  between  sorl  strength  properties  and  the  slip-shear  parame¬ 
ters  that  were  found  to  yield  good  simulation  of  the  experimental 
pull  coe L f ic ient-s lip  curves.  These  studies  showed  no  close  cor¬ 
relation  between  slip-shear  and  strength  parameters. 

Further  studies  indicated  that  the  relation  of  tire  dimensions 
to  load  has  some  effect  on  the  slip-shear  parameters.  A  correlation 
was  found  to  exist  between  the  dimensionless  numerics  developed  at 
WES  and  the  slip-shear  parameters.  The  following  relationships  were 
tentatively  established: 

Yuma  sand : 

j  =  0.046  -  0.0006  N 

J  o  s 

(7) 

-0.58  x  log  2N& 

K  =  e 


where 


N 

s 


CGR .  (b  .  d) 


3/2 


W 


h 


39 


IX.  PREDICTION  OF  TIRE  PERFORMANCE 
BY  THE  PROPOSED  TIRE -SOIL  MODEL 


In  Section  VIII,  relationships  were  established  that  allow 
the  prediction  of  tire  performance  from  input  parameters  that  are 
generally  available  or  can  be  obtained  by  conventional  testing 
techniques.  To  evaluate  the  accuracy  of  performance  predictions 
using  these  relationships,  performance  calculations  were  made  using 
the  computer  program  described  in  Appendix  B  for  all  cases  for 
which  experimental  data  were  available.  It  should  be  emphasized 
that  the  comparison  of  experimental  and  predicted  performance 
shown  graphically  in  the  subsequent  figures  represents  results 
that  can  be  obtained  from  data  generally  available  for  mobility 
eva lua  tion . 

Figures  10  through  36  show  predicted  and  measured  pull  coeffi¬ 
cients  as  a  function  of  slip  for  tests  1  through  27  perlormed  in 
Yuma  sand.  Figures  37  through  47  show  the  same  for  tests  Cl-1 
through  Cl-ll  performed  in  Buckshot  clay. 

Figures  48  through  30  show  typical  tire  centerline  geometries 
and  outlines  of  slip  line  fields  as  obtained  in  the  computer  pro¬ 
gram  for  the  various  slip  rates  in  tests  1,  5,  and  IS,  respectively. 
Figure  51  shows  the  same  in  the  case  of  a  strong  cohesive  soil  when 
there  is  only  one  rearward  failure  zone  (test  Cl-1). 

Figures  52  through  56  show  comparisons  of  predicted  and  mea¬ 
sured  torque  coefficients  in  some  typical  cases.  Note  that  the 
torque  coefficient  is  not  very  well  defined  for  the  centerline 
geometry  of  the  tire-soil  model  and  questions  arise  as  to  which 
value  of  the  continuously  changing  tire  radius  should  be  used  in 
the  denominator  of  the  formula  for  the  torque  coefficient.  To 
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v«y r«pv  y,  7CV-"  '''>J.' 


".  y^»T'«??>T»vi-w»|^'  .  ■ 


avoid  this  ambiguity  the  nominal  tire  radius  was 
putation  of  the  torque  coefficient  from  both  the 
computed  torque. 


used  in  the  com- 
mea  ured  and 


Figures  57  and  58  show  comparisons  of  measured  and  predicted 
slnkages  in  some  typical  cases. 
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Fig.  10  Pull  Coefficient  Versus  Slip  -  Test  1 
Yuma  Sand,  CGR  =  8.9 
Tire:  9.00-14,  Load:  440  lb 
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Pig.  11  Pull  Coefficient  Versus  Slip  -  Test  2 
Yuma  Sand,  CGR  =  16.4 
Tire:  6.00-16,  Load:  240  lb 
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Fig.  13  Pull  Coefficient  Versus  Slip  -  Test  4 
Yuma  Sand,  CGR  <=  15.6 
Tire:  6.00-16,  Load:  450  lb 
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1  ig .  14  Pull  Coefficient  Versus  Slip  -  Test  5 
Yuma  Sand,  CGR  -  13.0 
Tire:  4.00-7,  Load:  210  lb 
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Fig.  16  Full  Coefficient  Versus  Slip  -  Test  7 
Yuma  Sand,  CGR  «=  19.6 
Tire;  4.00-7,  Load:  220  lb 


0.8 


f 

I 


*  Predicted 
+  Measured 

I 


I 

I 


Pull 

Coefficient 


-20 


Slip,  % 


80 


Fig.  17  Pull  Coefficient  Versus  Slip  -  Test  8 
Yuma  Sand,  CGR  «=  19.9 
Tire:  4.00-7,  Load:  440  lb 
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Fig.  18  Pull  Coefficient  Versus  Slip  -  Test  9 
Yuma  Sand,  CGR  =  13.8 
Tire:  9.00-14,  Load:  840  lb 
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Fig.  20  Pull  Coefficient  Versus  Slip  -  Test  11 
Yuma  Sand,  CGR  =  12.4 
Tire:  9.00-14,  Load:  440  lb 
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Fig.  21  Pull  Coefficient  Versus  Slip  -  Test  12 
Yuma  Sane,  CGR  =22.5 
Tire:  9.00-14,  Load:  210  lb 
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Fig.  23  Pull  Coefficient  Versus  Slip  -  Test  14 
Yuma  Sand,  CGR  =  12.4 
Tire:  9.00-14,  Load:  870  lb 
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Fig.  24  Pull  Coefficienc  Versus  Slip  -  Tes 
Yuma  Sand,  CGR.  *=  13.0 
Tire:  9.00-14,  Load:  225  lb 
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’  Pull  Coefficient  Versus  Slip  -  Test  16 
Yuma  Sand ,  CGR  <=  24.2 
lire:  9.00-14,  Load:  670  lb 
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Fig.  26  Pull  Coefficient  Versus  Slip  -  Test  ln 
Yuma  Sand,  CGR  =  14.4 
Tire:  4.00-20,  Load:  225  lb 
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Fig.  27  Pull  Coefficient  Versus  Slip  -  Test  18 
Yuma  Sand,  CGR  =  12.0 
Tire:  4.00-20,  Load:  450  lb 
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Fig.  28  Pull  Coefficient  Versus  Slip  -  Test  ..9 
Yuma  Sand,  CGR  ■=  16.1 
Tire:  4.00-20,  Load:  330  lb 
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29  Pull  Coefficient  Versus  Slip 
Yuma  band,  CGR  =  16.4 
Tire:  4.00-20,  Load:  440  lb 
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I-ig.  it  Pull  Coefficient  Versus  Slip  -  Tes 
Yuma  Sand ,  CGR  =  3.2 
Tire:  4.00-20,  Load:  240  lb 
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Fig.  33  Pull  Coefficient  Versus  Slip  -  Test  24 
Yuma  Sand ,  CGR  =  13.0 
Tire:  9.00-14,  Load:  850  lb 
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Fig.  35  Puli  Coefficient  Versus  Slip  -  Test  26 
Yuma  Sand,  CGR  =  15.0 
Tire:  6.00-16,  Load:  220  lb 
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!•'  i  p .  38  Pull  Coefficient  Versus  Slip  -  Test  CL-2 
Buckshot.  Clay,  Cl  -  22 
Tire:  U.()0-J4,  Load:  h50  lb 
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» 1  Pu ] 1  doe  f  L . cienL  Versus  Slip  -  Test  C  L -  5 
Buckshot  Ciay,  Cl  «=  37 
Tire:  ft. 00 -lft,  Load:  8<S0  lb 
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ig.  44  Pull  Coefficient  Versus  Slip  -  Test  CL-8 
BucksnoL  Clay,  Cl  =  16 
Tire:  16  x  15-  \  Load:  450  lb 
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Fig.  52  Torque  Coefficient  Versus  Slip  -  Test  5 
Yuma  Sand,  CGR  «  13.0 
Tire:  4.00-7,  Load:  210  lb 
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Fig.  53  Torque 
Yu»t-  S£ 


Torqu</ Coefficient  Versus  Slip 
Yu^gc-  Sand,  CGR  °=  10.6 

re:  4.00-7,  Load:  220  lb 
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Fig.  54  Torque  Coefficient  Versus  Slip  -  Test  34 
Yuma  Sand,  CGR  «=  12.4 
Tire:  9.00-14,  Load:  870  lb 
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Fig.  55  Torque  Coefficient  Versus  Slip 
Yuma  Sand,  CGR  >=13.0 
Tire:  4.00-30,  Load:  450  lb 
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Fig.  56  Torque  Coefficient  Versus  Slip  -  Test  19 
Yuma  Sand,  CGX  =  16.1 
Tire:  7.00-20,  Load:  330  lb 
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iig.  57  Sinkage  Versus  Slip  -  Test  14 
Yuma  Sand,  CGR  =  12.4 
Tire:  9.00-14,  Load:  810  lb 
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Fig.  58  Sinkage  Versus  Slip  -  Test  17 
Yuma  Sand,  CGR  ■=  14.4 
Tire:  4.0U-20,  Load:  225  lb 
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X.  EVALUATION  OF  PERFORMANCE  PREDICTIONS 


The  comparison  of  computed  and  measured  performance  coeffi¬ 
cients  presented  in  Section  IX  shows  a  predictive  capability  of  the 
tire-soil  model  that  is  truly  remarkable  in  view  of  the  many  assump¬ 
tions  that  were  employed  to  make  the  model  simple.  It  was  expected 
that  three  dimensional  effects  and  bow  wave  effects  that  may  be  ap¬ 
preciable  under  certain  c ircums t mtes  and  were  not  accounted  for  in 
the  mode L  would  have  caused  greatei  tter  in  the  predictions  than 
actually  occurred.  These  unaccounted-for  effects  were  analyzed  in 
detail.  Results  of  the  analyses  are  presented  in  Section  XII. 

In  all  computations  a  uniform  angle  o,  interface  friction,  5, 
was  assumed.  This  assumption  resulted  in  inaccuracies  of  perfor¬ 
mance  predictions  for  the  rigid  wheel  (Ref.  13)  whenever  6  showed 
an  appreciable  variation  over  the  contact  area.  The  performance 
predictions  by  the  tire-soil  model  do  not  appear  to  be  affected  by 
possible  deviations  from  the  assumed  uniform  distribution  of  6. 

The  most  likely  explanation  for  this  is  that  the  distribution  of 
6  is  probably  close  to  a  uniform  one  in  the  case  of  tires.  The 
deformability  of  the  tires  and  their  ability  to  comply  to  the  kine¬ 
matic  conditions  at  the  interface  are  very  likely  to  result  in  a 
uniform  straining  and  associated  interface  friction  over  the  con¬ 
tact  area. 

An  interesting  feature  of  the  tire-soil  model  is  that  it  shows 
the  drop  in  the  pull-coefficient  with  slip  increasing  past  the  op¬ 
timum.  The  shape  of  the  pull-coefficient  versus  slip  curve  in 
frictional  soils  is  similar  to  that  of  the  shear  stress-displace¬ 
ment  curve  in  direct  shear  stress  tests  experienced  in  dense  soils. 
The  decrease  in  pull-coefficient  was  often  attributed  to  the  decrease 
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in  mobilized  shear  stress  with  lip.  However,  he  mobilized  shear- 
slip  curve  used  with  the  present,  l  n  e-soil  modi  L  and  rep  isenLto  b> 
Eq •  (5)  does  not  show  a  peak  but  increases  mono ionic a  1 1 y .  Thus  t no 
decrease  of  the  pu 1 1 -coef fie ieni  at  high  s  lip  rale  cannot  ne  at¬ 
tributed  to  this  effect,  buL  is  the  conseque-  of  tire-scd  1  :n~ 
teraction.  When  the  shear  stress  at  Lhe  interface  approaches  j ls 
maximum  value,  Lhe  normal  stresses  in  the  rear  field  thaL  an 
trolled  by  so1'.!  failure  conditions  rapidly  decrease.  The  tire- 
soil  model  responds  to  this  decrease  by  an  increase  of  the  entry 
and  rear  angles  and  of  the  sinkage  associated  with  these  angles. 

The  average  inclination  of  the  contact  area  increases  as  does  Lhe 
horizontal  component  of  the  normal  stresses  acting  on  this  area. 

The  net  result  is  a  decrease  of  the  pul L-coef f ic i ent  even  though 
the  mobilization  of  Lhe  shear  strength  at  the  interface  is  close 
to  its  maximum. 

The  above  finding,  is  significant  inasmuch  as  strain  softening 
of  the  interface  shear  strength  may  be  discounted  as  the  cause  of 
drawbar  pu 1 L  decrease  at  high  slip  rates.  Strain  softening  is 
characteristic  of  the  shear  strength  properties  of  dense  granular 
materials  where  part  of  the  shear  strength  is  derived  from  the 
necessity  that  particles  override  each  other.  When  Lhe  dense  mate¬ 
rial  loosens  up  shear  may  occur  without  much  override,  and  the  shear 
strength  decreases.  At  the  interface  the  mechanism  of  shear  strengtl 
mobilization  is  different  from  that  described  above  in  that  the 
failure  surface  is  a  solid  one  where  shear  may  occur  without  par¬ 
ticle  overriding.  Thus  the  monotonical  increase  of  shear  strength 
mobilization  as  intuitively  suggested  by  Janosi  and  Hanamoto  and 
represented  by  Eq.  (5)  is  more  representative  of  the  mechanism  of 
shear  at  the  interface  than  other  relationships  that  show  a  peak. 


While  the  decrease  in  pull-coefficient  was  well  reproduced  by 
the  model,  the  accompanying  significant  increase  in  the  torque  co- 
efrii  enc  was  not  always  reproduced.  Figures  52  through  54  show 
an  ir  rease  of  the  measured  torque  at  high  rates  of  slip  while  the 
predi  ed  torque  decreases.  On  the  other  hand,  Fig.  55  shows  a 
case  where  both  the  measured  and  predicted  torque  coefficient  de¬ 
crease  at  high  slip  rates.  An  analysis  ot  the  possible  causes  of 
this  discrepancy  was  made  and  the  followim,  com  lusions  were  reached. 

Th  .ire-soil  model  does  nut  account  r  the  internal  friction 
caused  by  the  deformation  of  arcass  and  fr  ctlon  in  the  bearings 
whi  le  the  effect  of  these  are  included  in  the  torque  measurements., 

It  is  estimated  that  discrepancies  on  this  account  ma\  unount  to 
2-3  percent  in  the  value  of  the  torque  coeLt  orient  . 

Some  rotational  inertia  <  fleets  may  have  affect  id  t  he  Lorque 
measurements.  Inaccuracies  on  this  account  are  probably  small. 

The  centerline  geometries  and  radius  variations  assumed  in 
the  tire-soil  model  for  these  high  slip  rates  were  .  estivated. 

It  was  found  that  differences  in  Lorcin  prediction  on  acco  nt  of 
misestimated  radii  are  not  likely  to  result  in  differon  e  >f  more 
than  2-3  percent  in  t  ie  predicted  value  of  the  ti  rque  coeiiicient. 

The  assumed  smooth  distribution  of  interlace  stresses  in  con¬ 
trast  to  that  observed  in  experiments  that  shoe  peaks  due  to  car¬ 
cass  stiffness  may  result  in  underestimated  Lorque.  The  peaks  oc¬ 
cur  close  Lc  the  entry  and  exit  angles  where  deflection  is  less  and 
radius  is  greater  than  in  the  center  (Fig.  59).  Again,  differences 
on  this  account  are  not  likely  to  exceed  2-3  percent  in  the 
torque  value  . 

Even  if  the  sources  of  discrepancies  act  cumulatively,  the 
over-all  effect  wourd  not  be  more  than  10  percent  in  the  value 
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of  torque  coefficient.  Since  at  very  high  slips  discrepancies  oc¬ 
curred  in  the  torque  coefficient  as  high  as  20  percent,  there 
must  be  yet  another  source  of  these  discrepancies.  It  was  found 
that  the  most  likely  cause  of  discrepancies  higher  than  10  per¬ 
cent  is  that  the  assumption  that  the  moment  of  interface  normal 
stresses  about  the  axle  is  zero  does  nr  t  hole*  true  at  high  slip 
rates.  This  assumption  was  based  on  the  experimental  finding  of 
Fret  tag  et  al.  (Ref.  4).  It  appears  that  generalization  of  the 
validity  of  the  conclusion  drawn  in  Ref.  4  from  a  limited  number 
of  experiments  (that  the  moment  of  normal  stresses  is  zero)  wa^  not 
warranted.  It  appears  that  at  high  slip  rates  the  center  of  tire 
axis  shifts  in  such  a  way  that  normal  interface  stresses  gene1  ate 
moments  about  the  axis.  Such  a  shift  in  the  relative  position  of 
the  axis  to  the  interface  was  observed  in  the  experiments  -where 
centerline  geometry  was  measured  (Fig.  60).  Unfortunately,  the^e 
were  no  simultaneous  interface  stress  measurements  in  these  ex¬ 
periments  so  that  there  is  no  way  1 n  check  the  validity  of  the 
assumption  that  the  torque  from  normal  stresses  is  zero. 

The  specif  ative  conclusion  that  at  high  slip  rates  interface 
normal  stresses  may  cause  moments  about  the  axis  is  indirectly 
supported  by  the  dimensional  analysis  of  tire-soil  interaction 
(Ref.  13)  and  the  comparatively  poor  torque  p  rformance  prediction 
with  the  WES  method.  In  the  dimensional  analysis  it  was  found  that 
the  relation  of  the  torque  number  to  the  sand  number  was  not  well 
defined  (p.  11/ ,  Ref.  18)  and  torque  performance  predictions  were 
found  to  separate  for  various  tire  shapes  (p.  15,  Ref.  18).  Ob¬ 
viously,  a  shift  in  the  relative  position  of  the  Lire  axis  to  the 
interface  is  a  tire  property  that  cannot  possibly  be  cl  acterized 
by  a  single  deflection  number  obtained  under  vertical  loading  con¬ 
ditions  on  a  rigid  surface  and,  therefore,  cannot  be  accounted  for 
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by  dimensional  analysis  techniques  using  this  tire  characceriza- 
tion.  A  shift  in  the  position  of  the  axis  would  affect  the  torque 
coefficient  hut  not  the  other  performance  cha rac Leris  tics  and 
would  explain  the  discrepancies  of  torque  coefficient  predictions 
by  both  the  tire-soil  model  and  WES  techniques. 

The  performance  prediction  for  Lest  13  was  analyzed  in  detail 
since,  in  contrast  to  Lhe  other  tests,  performance  prediction  in 
this  case  was  relatively  poor  (Fig.  22) .  A  look  at  the  pull  per¬ 
formance  prediction  curve  established  by  WES  (Fig.  61)  reveals 
that  this  curve  is  very  steep  at  the  sand  numeric  for  test  13 
(N^  =  5.06).  This  indicates  that  under  these  conditions  a  small 
variation  in  soil  strength  results  in  a  large  variation  of  the 
pull-coefficie.it . 

An  investigation  was  made  to  determine  the  effect  of  such 
small  variations  on  the  tire  performance  predicted  by  the  model. 
Figure  22  shows  pull-coefficients  computed  for  i  =  37.6°,  es¬ 
tablished  by  Eq .  (3),  and  for  ■:  =  36°  and  35°.  These  two  lat¬ 
ter  friction  angles  yield  acceptable  performance  prediction  and 
are  within  the  limits  of  accuracy  by  which  the  friction  angle  can 
be  determined  from  cone  penetration  measurements.  Thus  modifica¬ 
tion  of  the  tire-scil  model  on  account  of  Lest  13  was  not  deemed 
necessary. 

The  pu 1 1 -coe f J io ien t  prediction  for  test  Cl-1  was  also  analyzed 
in  detail.  These  analyses  indicated  that  it  was  not  possible  to 
duplicate  the  pu 1 1 -coe f f ic ien t  of  0.56  observed  in  the  test  at 
60  percent  slip  with  any  reasonable  variation  of  the  model  parame¬ 
ters.  An  examination  of  the  predicted  sinkages  in  these  analyses 
indicated  that  at  20  percent  slip  (where  sinkage  measurement  was 
available)  the  predicted  sinkage  was  less  than  the  measured  one 
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excluding  the  possibility  that  underestimation  of  the  pull  coeffi¬ 
cient  resulted  from  an  overes tima ted  sinkage.  For  small  sinkages 
the  contact  area  may  be  approxina ted  by  a  plane  inclined  at  an 
angle  e  to  the  horizontal,  and  the  following  approximate  rela¬ 
tionships  for  the  maximum  value  cf  the  pull-coefficient  may  be 
developed  (Fig.  62). 

T  =  N  •  tan  r  +  c  .A  (9) 

max 

H  =  (T  -  T  )  cos  e  =  T  •  cos  c  -  W  tan  c  (10) 
max  max  max 


A  =  N/p 


av 


(il) 


Pull-coefficient  =  (H/W)  =  .tan  ;  + 

rna  x  \ 


av 


cos  c  -  tan 


(12) 


An  examination  of  Eq.  (12)  shows  that  for  the  pull-coefficient 
to  become  as  high  as  0.56,  either  tan  ;  or  c/p  or  both 
would  have  to  have  much  higher  values  than  assumed  or  calculated 
in  the  model.  It  is  highly  unlikely  that  the  clay  as  prepared 
for  the  test  at  a  high  degree  of  saturation  would  exhibit  a  fric¬ 
tion  angle  appreciably  higher  than  13° „  Ln  cohesive  soils  the 
•/alue  of  c/p  is  not  appreciably  affected  by  the  value  of  co¬ 
hesion  when  the  soil  strength  controls  the  interface  stresses. 

In  such  a  case  p  is  approximately  five  times  the  cohesion  re¬ 
sulting  in  a  c/p  value  of  about  0.2.  The  value  of  c/p 

could  be  significantly  higher  only  if  the  limiting  pressure,  p^, 
controlled  the  average  pressure  p  .  In  order  for  this  Lo  be  the 
case,  the  actual  inflation  pressure  would  have  had  to  be  much  lower 
than  the  design  pressure.  Although  it  is  conceivable  that  some 
leakage  occurred  during  the  test,  it  is  unlikely  that  a  signifi¬ 
cant  drop  in  the  inflation  pressure  would  have  gone  unnoticed. 
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A  factor  Lhal  may  increase  Lhe  pu 1 i -c oe i f ic ien  L  and  is  not 
accounted  for  in  the  model  is  the  effect  of  side  wa 1 L  adhesion. 

An  estimate  of  the  effect  of  this  adhesion  showed  that  it  could 
not  affect  the  pul  1 -coef f ic ien t  by  more  than  2-3  percent. 

All  the  above  considerations  indicate  that  the  underestimation 
of  the  pull-coefficient  ^.n  this  case  cannot  be  attributed  to  an  un¬ 
satisfactory  simulation  by  the  model  or  to  inaccuracies  in  the 
estimates  of  the  strength  properties  of  soil.  The  most  likely  ex¬ 
planation  for  the  discrepancy  between  observation  and  model  simula¬ 
tion  is  that  soil  conditions  were  not  ui.’ form  during  Lhe  test  and 
the  strength  of  soil  at  the  surface  was  higher  than  the  strength 
averaged  over  a  six-inch  depth.  In  experiments  with  clays  close 
to  saturation  it  is  very  difficult  to  prevent  evaporation  from  Lhe 
surface;  a  relative  humidity  corresponding  to  Lhe  moisture  equi¬ 
librium  of  the  soil  would  have  to  he  maintained  during  the  prepa¬ 
ration  of  the  soil  bed  as  well  as  during  the  Lest.  In  earlic'  ex¬ 
periments  at  WES  a  decrease  in  the  moisture  content  of  the  upper¬ 
most  layer  in  the  soil  bed  was  indeed  observed  (Ref.  19),  and  at¬ 
tention  was  called  to  the  importance  surface  strength  in  Lire 
performance  tests.  Unfortunately,  <  penetrometer  tests  are  noi. 
very  well  suited  to  discriminate  between  surface  soil  strength  and 
the  average  strength  of  a  six-inch  deep  layer,  and  later  investiga¬ 
tions  were  confined  to  Lhe  detrimental  effect  of  weak  (wet)  surface 
layers  . 

The  strength  of  Lhe  surface  layer  affects  the  pu 1 1 -coe f f 1 cient 
in  the  following  way:  In  Eq.  (9)  for  the  maximum  traction,  Lh 
strength  parameters  of  the  surface  layer  may  be  used.  The  failure 
conditions  in  the  sol,  however,  are  controlled  by  the  lower  strength 
of  the  soil  below  the  surface  layer  so  that  the  normal  pressures  and 


p  remain  essentially  unaffected  by  the  higher  strength  of  the 
surface  layer.  Thus  in  Eq.  (12)  both  m  and  c  refer  to  the 
stronger  surface  layer  while  p  corresponds  to  the  lower  strength 
of  the  underlying  layer.  The  result  is  a  significantly  higher  pull- 
coefficient  than  that  obtainable  unde1'  uniform  soil  conditions. 

The  potential  effect  of  a  strong  surface  layer  on  the  pull- 
coefficient  is  subject  to  limitations  imposed  by  the  underlying 
layer.  The  strength  of  the  upper  layer  to  develop  traction  can 
only  be  utilized  to  the  extent  that  the  shear  stresses  at  the  in¬ 
terface  between  the  upper  and  lower  la^er  do  not  exceed  the  failure 
criteria  for  the  lower  layer.  This  will  depend  on  the  thickness  of 
the  upper  layer  and  the  relative  strength  of  the  upper  layer  tc 
that  of  the  lower  one.  These  conditions  may  be  analyzed  by  ex¬ 
panding  the  present  tire-soil  model  so  as  to  allow  for  a  two-layer 
system  or  for  soil  strength  changing  gradually  with  depth.  Such  r 
model  would  be  essentially  the  same  as  proposed  previously  for  the 
analysis  of  slipperiness.  Slipperiness  is  essentially  the  same 
phenomenon  as  discussed  previously  except  that  the  strength  ratios 
of  the  upper  and  lower  layer  are  reversed. 


XI.  SIMULATION  OF  TIRE  PERFORMANCE  IN  COHESIVE -FRICTIONAL  SOILS 


The  comparisons  of  predicted  tire  performances  with  measured 
ones  given  in  Section  X  show  that  the  proposed  tire-soil  model 
yields  very  good  predictions  for  a  variety  of  tires  and  soil  con¬ 
ditions.  The  experimental  data  available  for  these  comparisons 
however,  were  restricted  to  two  types  of  soil:  the  purely  fric¬ 
tional  Yuma  sand  or  the  purely  cohesive  Buckshot  clay.  From  the 
viewpoint  of  the  tire-soil  model  this  restriction  is  not  considered 
as  a  limitation  to  the  use  of  the  model  in  frictional-cohesive 
soils.  On  the  contrary,  the  computation  of  interface  stresses  by 
plasticity  theory  takes  into  consideration  the  simultaneous  effect 
of  friction  and  cohesion  and,  therefore,  it  can  be  assumed  that  the 
tire-soil  model  simulates  tire  performance  in  cohesive -fric t iona I 
soils  as  well  as  in  the  cases  shown  in  Section  X.  Nevertheless, 
it  seemed  desirable  to  evaluate  model  performance  under  soil  con¬ 
ditions  other  thin  tested  at  WES.  The  tire  performance  tests  con¬ 
ducted  at  the  National  Tillage  Machinery  Laboratory  in  various 
types  of  agricultural  soil  appeared  to  be  suitable  for  this  pur¬ 
pose.  However ,  these  tests  were  performed  for  comparative  tire 
performance  studio?  where  soil  conditions  were  intended  to  dupli¬ 
cate  conditions  likely  to  be  encountered  in  the  field  by  agricul¬ 
tural  tractors  (for  example,  moisture  content  increasing  with  depth). 
Since  at  this  time  the  tire-soil  model  has  not  yet  been  expanded  to 
include  layered  or  nonuniform  soil  conditions,  the  simulation  of 
the«e  agricultural  tire  performance  tests  by  the  model  strictly 
applicable  to  homogeneous  soils  can  only  be  considered  as  a  crude 
approximation . 

The  tire  performance  tests  selected  for  simulation  were  part 
of  a  larger  program  where  the  effect  of  lug  angle  on  performance 
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was  investigated.  These  tests  are  described  in  detail  in  Ref.  20. 

The  selected  tests  are  the  ones  performed  with  smooth  tires  that 
were  tested  in  the  large  program  for  reference  purposes. 

In  order  to  have  performance  evaluations  comparable  to  that 
reported  in  the  preceding  sections,  the  digital  data  points  re¬ 
ceived  on  a  punched  paper  tape  from  the  National  Tillage  Machinery 
Laboratory  for  various  travel  reduction  rates  had  to  be  converted 
to  slip  rates.  In  tnis  conversion  in  the  absence  of  tire  deflection 
data,  the  so-called  loaded  tire  radius  taken  from  the  manufacturer's 
catalog  was  used  in  the  computation  of  the  slip  rate.  The  slip 
rates  shown  in  the  various  figures  may  be  in  a  slight  error  on  this 
account.  For  each  test  over  100,  data  points  were  obtained. 

Figure  63  shows  data  points  for  test  A1  Logether  with  the  best 
fitting  second  degree  polynomial  obtained  by  the  appropriate  com¬ 
puter  library  program  and  displayed  on  the  visual  display  terminal. 

In  the  following  comparisons  of  performance  predictions  with  National 
Tillage  Machinery  Laboratory  test  results  only  the  best  fitting 
curves  are  shown  to  preserve  claiity.  The  best  fitting  curves  for 
tests  Al  and  A2  shown  in  Fig.  64  and  those  for  tests  A3  through  A6 
shown  in  Fig.  65  are  almost  identical,  indicating  the  good  repro¬ 
ducibility  of  these  tests. 

Results  of  the  simulation  of  these  tests  by  the  tire-soil  model 
are  shown  in  Figs.  64  and  65  by  crosses.  The  results  of  the  tests 
performed  in  the  Decatur  silty  1 oam  were  duplicated  fairly  closely 
using  the  strength  properties  determined  by  triaxial  tests  for  the 
uppermost  layer  with  the  least  moisture  content.  Simulation  of 
the  tests  performed  in  the  Vaiaen  clay  is  not  as  good  as  in  the 
Decatur  silty  loam,  although  in  the  simulation  sham,  a  cohesion 
value  about  30  percent  higher  was  used  than  that  determined  by 
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triaxial  lests  for  the*  uppermost  layer  of  the  Vaideu  c  lay.  In 
tin  simulation  s Lip  parameters  were  assumed  by  trial  and  error 
proci  lure  . 

Il  ippears  LhaL  in  order  to  obtain  good  s  i  inn  1  a.  L  i  on  ,  soil 
strength  values  equal  or  higher  than  that  of  the  strongest  upper¬ 
most  layer  would  have  to  be  used  with  the  tire-soil  model.  One 
explanation  for  this  is  that  Lhe  soil,  strength  determination  by 
triaxial  Lests  v<,us  onducled  at  a  lower  rate  than  Lhe  loading  rale 
applied  to  the  soil  in  the  Lire  test.  Another  possibility  that  at 
first  glance  may  seem  paradoxical  is  that  in  a  layered  system  with 
strength  properties  decreasing  with  depth  a  higher  pull-coefficient: 
can  be  obtained  than  in  a  uniform  soil  that  would  exhibit  the 
sirength  properties  of  the  uppermost,  strongest  layer.  There  are 
indications  that  a  weaker  underlying  layer  would  act  in  a  similar 
way  as  a  low  inflation  pressure:  it  would  reduce  the  normal  in¬ 
terface  stresses  while  allowing  interface  shear  stresses  as  high 
as  the  strength  of  the  upper  layer.  The  extension  of  Lhe  tire- 
soil  model  to  layered  soil  system  and  nonuniform  soil  conditions 
would  allow  the  analysis  of  such  cases;  until  such  a  model  is 
available  no  conclusive  evaluation  of  the  test  results  can  be  made. 
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XII.  EVALUATION  OF  THE  TIRE-SOIL  MODEL  AND  SUGGESTED  IMPROVEMENTS 


The  prediction  of  tire  performance  presented  in  the  various 
figures  of  Section  IX  is  very  good  and  one  could  say  much  better 
than  might  be  expected  from  a  two  dimensional  model  that  cannot 
take  into  account  three  dimensional  effects.  Since  three  dimen¬ 
sional  effects  as  well  at  bow  waves  certainly  affected  tire  per¬ 
formance  at  least  in  some  of  the  tests  conducted,  it  was  '"hought 
necessary  to  analyze  these  effects  and  determine  why  it  was  pos¬ 
sible  to  achieve  good  prediction  accuracy  without  the  considera¬ 
tion  of  these  effects.  Following  are  .he  results  of  these  analyses. 

It  was  pointed  out  in  the  previous  report  01  rigid  wheel-sriJ 
interaction  (Ref.  13)  that  lateral  failure  conditions  constitute  a 
limitation  of  the  normal  stresses  that  may  arise  at  the  interface, 
The  interface  stress  measurements  performed  with  rigid  wheels  ei„ 
clusively  iridicateu  that  the  di..triDution  of  normal  stresses  across 
the  interface  is  no  longer  ur.i  orm  if  lateraL  failure  conditions 
prevail.  The  limitation  on  the  normal  stresses  imposed  by  lateral 
failure  conditions  may  be  formulated  and  computed  by  plasticity 
theory;  this  development  was  outside  the  scope  of  the  present  pro¬ 
gram.  However,  for  the  purpose  of  this  discussion  an  approximate 
formula  was  developed  for  the  computation  of  normal  stresses  at 
the  outside  edge  of  the  tire  on  the  basis  of  lateral  failure.  In 
the  case  of  cohesionless  soils  the  ma <imum  normal  stress  that  can 
arise  at  the  edge  of  the  interface  depends  on  the  friction  angle 
and  on  the  depth  of  the  point  of  the  interface  below  the  original 
surface.  Typical  edge  stresses  computed  by  this  approximate  for¬ 
mula  are  compared  with  normal  stresses  computed  in  the  tire-soil 
model  in  Figs.  66  through  68.  Centerline  geometries  for  the  same 
test  are  shown  in  Fig.  49.  It  can  be  seen  that  the  smaller  the 
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Fig.  68  Limitation  of  Interface  Normal  Stresses  Imposed  by  Lateral 
Failure  Conditions  -  Test  5,  36%  Slip,  Sinkage  =  1.9  in. 


sinkage,  the  more  pronounced  the  effect  of  lateral  failure  condi¬ 
tions.  In  the  case  of  small  sinkage,  the  pull-coefficient  is  not 
influenced  by  the  average  inclination  of  the  contacL  area  and  is 
essentially  determined  by  the  a / c  ratio.  In  such  a  case  the 

ri 

tire-soil  model  predicts  the  pul  1 -coefficient  accurately  even 
though  the  computed  arc  lengLh  of  the  front  field  is  snorter  than 
it  would  be  had  the  lateral  effect  been  considered.  With  increasing 
sinkage,  the  limitation  on  the  normal  stresses  due  to  lateral 
failure  becomes  less  significant  and  preuxction  is  acceptable.  In 
cohesive  soils,  the  effect  of  lateral  failure  is  less  significant. 
Even  though  this  fortuitous  situation  allows  a  .satisfactory  per¬ 
formance  prediction  by  the  two  dimensional  model,  consideration  of 
the  three  dimensional  effect  would  result  in  the  following  advantages 

®  In  the  present  two  dimensional  model  Lhe  re  la t i onshins 
established  for  the  estimate  of  e  are  differ  nL  for 
the  various  tire  types.  It  is  believed  that  a  similar 
model,  expanded  to  take  the  limitations  on  the  inter¬ 
face  normal  stresses  due  to  lateral  failure  into  con¬ 
sideration,  would  probably  result  in  a  uniform  relation 
for  e  for  all  tire  sizes. 

Consideration  of  lateral  failure  in  Lhe  model  would 
allow  performance  simulation  of  dual  and  multiple  tire 
arrangements  with  proper  consideration  of  their  spacing. 

0  The  development  of  a  technique  to  consider  lateral 
failure  would  be  directly  applicable  to  tracks  and 
track  elements  and  could  be  the  basis  of  tracked 
vehicle  modeling. 

As  regards  the  bow  wave  effect,  the  following  considerations 
apply:  The  tire-soil  interaction,  as  was  shown  previously  can  be 
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satisfactorily  simulated  by  the  model  in  which  the  soil  is  repre¬ 
sented  by  its  strength  properties.  Sinkage  as  well  as  stresses 
in  the  soil  beneath  the  tire  are  determined  without  consideration 
of  the  volumetric  changes  in  the  soil.  However,  the  total  volume 
change  of  the  soil  due  to  the  stress  system  generated  by  the  tire 
load  should  be  in  balance  with  the  total  volume  change  represented 
by  the  size  of  the  rut.  If  there  is  an  unbalance,  bow  waves  form. 

In  the  experiments  performed  at  WES  with  densely  packed  wooden 
rollers  (a  two  dimensional  soil  model),  tire-soil  interaction  re¬ 
quired  a  certain  sinkage  and  volume  change  that  could  not  be 
balanced  by  further  densi  f  ica  t  ion  of  the  wooden  rollers.  Conse¬ 
quently,  sizable  bow  waves  developed  in  these  tests  even  at  high 
slip  rates.  Since  the  volume  change  of  the  soil  model  was  prac¬ 
tically  nil,  and  no  lateral  displacement  occurred  in  the  two  dimen¬ 
sional  model,  no  volume  balance  could  be  achieved  and  the  size  of 
the  bow  wave  did  not  stabilize  during  Lhe  test.  In  actuality  the 
conditions  are  never  rigorously  two  dimensional  and  significant 
bow  waves  occur  only  with  towed  and  relatively  wide  tires. 

The  effect  of  bow  waves  on  tire-soil  interaction  may  be  ana¬ 
lyzed  by  the  model  by  assuming  that  the  geometry  of  the  free  sur¬ 
face  in  the  front  of  the  tire  is  that  of  the  bow  wave.  Preliminary 
studies  indicated  that  it  is  feasible  to  account  for  the  bow  wave 
geometry  in  the  computation  of  the  front  slip  line  field.  Figure  69 
shows  the  slip  line  field  geometry  for  an  assumed  bow  wave.  Pre¬ 
liminary  analyses  indicate  that  the  effect  of  a  bow  wave  on  inter¬ 
face  stresses  and  tire  performance  is  not  as  great  as  the  physical 
appearance  of  the  bow  wave  would  suggest. 

While  it  would  be  desirable  to  take  the  effect  of  bow  waves 
into  account  in  tire-soil  interaction,  this  would  necessarily  in¬ 
volve  the  determination  of  at  least  another  parameter  that  expresses 
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Lhe  volume  change  properties  of  soil.  (This  is  discussed  in  more 
detail  in  Appendix  /..)  Fortunately,  as  our  preliminary  analyses 
indicate,  the  how  wave  effect  does  not  influence  significantly 
the  pu 1 1 -coe f f ic ien t  determination  for  driven  tires.  However, 
for  towed  and  very  wide  tires  the  bow  wave  effect  may  not  be 
negligible.  The  tire-soil  model  may  be  expanded  Lo  include  Lhe 
consideration  of  bow  wave  effects,  but  further  research  is  needed 
to  develop  methods  for  the  estimate  of  soil  volume  changes,  volume 
change  balance,  and  bow  wave  size. 
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XIII.  USE  OF  THE  TIRE -SOIL  MODEL  AS  AN  ANALYTICAL  TOOL 


The  simulation  of  suit  soil  tire  performance  by  mathematical 
models  is  advantageous  not  only  for  performance  predictions  as 
discussed  in  Sections  IX  through  XII,  but  also  for  the  analysis 
and  understanding  of  tire-soil  interaction  phenomena.  In  Section  X 
some  examples  of  the  use  of  the  tire-soil  model  for  the  analysis  of 
tire-soil  interaction  phenomena  have  already  been  presented.  For 
example,  it  was  shown  that  a  drop  in  the  pull-coefficient  experienced 
at  high  slip  rates  occurs  in  frictional  soils  even  if  the  developed 
shear  stress  increases  monotonicaliy  with  slip.  Another  analysis 
in  Section  X  showed  the  significance  of  surface  layer  strength. 

The  tire-soil  model  is  preeminently  suited  for  the  analysis 
of  the  effect  of  changes  in  any  one  input  variable  on  tire  per¬ 
formance.  For  example,  the  maximum  drawbar  pull  that  a  tire  can 
exert  in  sand  can  be  analyzed  and  the  maximum  drawbar  pull  at 
various  tire  loads  determined.  Figure  70  shows  the  results  of  such 
an  analysis  for  a  9.00-14  tire  inflated  to  16.4  psi.  The  soil  is 
Yuma  sand  with  a  cone  index  gradient  of  13.8.  It  can  be  seen  that 
there  is  a  maximum  drawbar  pull  that  cannot  be  exceeded  by  increas¬ 
ing  the  tire  load  even  if  the  available  torque  is  unlimited. 

Another  way  of  presenting  the  effect  of  tire  load  on  pull 
performance  is  shown  in  Fig.  71,  which  shows  the  decrease  of  the 
pull-coefficient  as  the  tire  load  increases. 

Parametric  analysis  of  design  variables  can  be  easily  performed 
by  the  mathematical  tire-soil  model.  In  Fig.  72  the  pull-coefficients 
for  various  tire  diameters  are  shown,  all  other  variables  being  con¬ 
stant.  The  fact  that  larger  diameter  tires  work  better  in  sand  than 
smaller  ones  is,  of  course,  well  known.  The  mathematical  model  con¬ 
firms  this  fact  in  quantitative  terms  over  the  whole  range  of  slip. 
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T.  io  n  |>fc  i  i  '  v  Lorque  and  sinkagc  values  arc  also  availab, 
thf  computer  program.  Performing  t  ho  performance  computation: 
for  soil  ondi lions  representative  of  a  mission,  Lhe  design  can 
be  optimized  and  trade-oft  studies  can  be  made  for  military  de¬ 
ls  i  c  ns  . 

The  effect  of  infla  ion  pressure  on  pull  per  I o  nuance  is  shown 
in  Fig.  7'i .  The  in;  lotion  press  a"*-:,  correspond  to  Is,  ?  i,  and 
J  j  percent  deflection  (in  WPS  terminology),  re  pecLively.  fin 
analysis  confirms  the  well-knc.wn  beneficial  effect  of  lowering  the 
inflation  pressure  on  traction. 

The  effect  of  soil  properties  on  tire  performance  may  also  be 
easily  analyzed.  An  example  of  such  an  analysis  is  shown  in  Fig.  7 1 
where  pull  performances  lor  the  conditions  of  test  9  are  compared 
for  various  values  of  the  unit  weight  . ;  all  other  -onditions  are 
unchanged.  It  Is  seen  that  a  slight  increase  in  improves  tire 

performance  by  some*  three  percent  in  the  pu  1  1  -cor  .  f  i  c  ient  .  Such  a 
difference  may  remain  obscured  in  an  experimental  program  where 
da  *  a  scat’er  occurs  for  many  reasons  and  a  detailed  s  tat  i  L  i  cm  1 
analysis  would  be  necessary  to  show  Lhat  such  observed  difference^ 
are  statistically  significant.  In  an  analysis  by  a  mathematical 
model  ail  other  variables  (including  friction  angle)  are  exactly 
the  same,  and  the  obtained  results  clearly  show  the  influence  of 
changes  in  Lhe  selected  variable,  the  unit  weight.  The  effect  of 
unit  weight  is  more  pronounced  under  other  circumstances.  In 
Fig.  74  computed  pu 1 1 -coe f f icien ts  are  also  shown  for  the  case  of 
=  h 5  Ib/cu  ft  that  corresponds  to  the  submerged  unit  weight  of 
soil.  It  is  seen  that  the  submergence  affects  tire  performance 
appreciably  even  if  only  its  effect  on  unit  weight  is  considered. 

In  this  analysis  pore  water  pressures  were  assumed  to  be  neutral; 
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the  effect  ol  pore  water  pressures  that,  depending  on  the  density 
of  soil,  decreases  or  increases  the  strength  of  soil  may  override 
the  effect  of  unit  weight. 

Although  the  comparison  of  predicted  and  measured  puJ 1  per- 
ormance  presented  in  Section  IX  includes  a  wide  range  of  soil 
strength  it  is  of  intei  st  to  show  the  variation  of  pull  per¬ 
formance  when  only  the  strength  is  varied,  al  other  conditions 
being  equal.  Such  a  comparison  is  shown  in  Fig.  75.  Soil  strength 
as  1  i_g .  75  and  the  other  performance  analyses  indicate,  is  the  most 
crucial  single  factor  in  vehicle  mobility. 

Figure  76  shows  uhe  effect  of  a  small  cohesion  on  the  pull 
performance  curve.  A  cohesion  of  10  lb/sq  ft,  which  is  indeed 
very  small  and  could  be  the  result  of  a  slight  dampness  in  fric¬ 
tional  soils,  improves  the  pull  performance  appreciably.  The  ef¬ 
fect  of  cohesion  on  pull  performance,  if  it  acts  in  conjunction 
with  a  high  i fiction  angle,  is  great. 

As  these  examples  indicate,  the  model  is  veil  suited  for  para¬ 
metric  analysis  of  various  design  features.  These  analyses  can  be 
performed  at  a  fraction  of  the  cost  of  experiments  and  in  a  matter 
of  minutes.  The  availability  of  such  a  model  is  expected  to  spur 
interest  in  optimizing  the  design  of  off-road  vehicles. 
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Fig.  70  Maximum  Drawbar  Pull  at  Various  Tire  Loads. 
Tire:  9.00-14,  Infl  Pressure:  16  4  psi 
Yuma  Sand,  CGR  =■•  13.8 
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Fig.  73  Effect  of  Inflation  Pressure  on  Pull  Performance 
Yuma  Sand,  CGR  «=  13.8 
Tire:  9.00-14,  Load:  830  lb 
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Fig.  74 


Effeet  of  Soil  Unit  Weight  on  Pull  Performer 
Yuma  Sand,  CGR  -  13.8,  Load:  830  lb 
Tire:  9.00-14,  Infl.  Pressure:  16.4  psi 
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Fig.  75  Effect  of  Soil  Strength  on  Pull  Performan< 
Yuma  Sand,  Load:  830  lb 
Tire:  9.00-14,  Infl.  Pressure:  16.4  psi 
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Fig.  7b  Effect  of  a  Small  Cohesion  on  Pull  Performance 
in  Sand,  Yuma  Sand,  Load:  830  lb 
Tire:  9.00-14,  Infl.  Pressure:  lb. 4  psi 
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XIV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  ma  Lhema  lica  1  Lirc-soil  inode  J  has  been  developed  Lhal  re¬ 
flects  the  role  of  Lire  deflection  and  soil  s  rength  in  tire-soil 
interaction  with  sufficient  fidelity  so  that  prediction  of  Lire 
performance  by  the  mathematical  model  is  within  the  range  of  the 
accuracy  of  experiments.  On  the  basis  of  the  mathematical  model 
a  computer  program  was  prepared  Lot  the  prediction  and  analysis 
of  tire  performance.  The  computer  program  requires  a  9K  core 
on  a  PDP  10  computer;  average  running  time  is  30-60  seconds. 

The  computer  program  can  also  be  accommodated  to  minicomputers; 
considerable  savings  in  running  cost  can  be  thereby  achieved. 

Model  input  parameters  are  of  three  kinds:  tire,  soil,  and 
slip  parameters.  For  the  determination  of  these  parameters  from 
conventionally  obtained  data  theoretical  studies  and  correlation 
analyses  have  been  made  with  the  following  results: 

6  For  the  tire  input  parameter,  p^,  a  generally  valid 
relationship  with  p^ ,  the  inflation  pressure,  has 
been  established.  The  input  tire  deflection  parameter, 
r  ,  has  been  related  to  the  deflection  on  rigid  surface 
by  various  relationships,  each  valid  for  a  particular 
Lire  type.  For  the  establishment  of  a  generally  valid 
relationship,  three  dimensional  effects  would  have  to 
be  considered  in  the  model.  It  is  recommended  that  the 
tire-soil  model  be  expanded  to  account  for  this  three 
dimensional  effect.  Such  an  expansion  would  also  allow 
the  meaningful  analysis  of  performance  of  dual  and 
multiple  tire  arrangements  with  respect  to  tire  spacing. 
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For  the  laboratory  determination  of  soil  input  parame¬ 
ters  a  triaxial  testing  technique  has  been  tentatively 
developed  that  varies  the  chamber  pressure  with  strain 
and  simulates  thereby  the  strain  and  stress  paths  in 
the  soil  under  tire  loads  better  Lhan  conventional  tech¬ 
niques.  It  is  recommended  that  further  research  be  con¬ 
ducted  for  the  establishment  of  strain  and  stre  path 
criteria  in  triaxial  testing  for  mobility  purposes  and 
appropriate  chamber  pressure  control  devices  he  developed. 


0  For  the  determination  of  soil  input  parameters  from  cone 
penetration  data,  appropriate  relationships  have  been  de¬ 
veloped  and  incorporated  with  the  computer  program.  These 
relationships  are  valid  for  either  frictional  or  cohesive 
soils.  For  f ric t iona i-cohes i ve  soils  it  is  recommended 
that  evaluation  methods  for  the  determination  of  soil 
strength  properties  from  field  tests  he  developed,  the 
sui Labi  lily  of  the  various  techniques  for  soil  strength 
determination  he  investigated,  and  appropriate  modifica¬ 
tions  in  these  techniques  undertaken. 

For  the  input  slip  parameters  correlations  with  the  sand 
or  clay  numeric  have  been  established  for  the  Yuma  sand 
and  Buckshot  clay,  respectively,  the  two  soils  for  which 
Lire  performance  data  were  available.  It  is  recommended 
that  a  theoretical  research  on  the  interface  shear  stress 
development  with  slip  be  initiated,  and  methods  for  the 
determination  of  these  two  parameters  in  ihe  laboratory 
and  in  the  field  be  developed. 

The  analysis  of  tire  performance  tests  in  clay  showed  that 
surface  conditions  may  affect  Lhe  pull  coefficient  appreciably,  not 
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only  if  the  surface  is  wetter,  but  also  when  it  is  drier  than  the 
soil  beneath  the  surface  layer.  Wet  surface  conditions  have  been 
of  concern  because  of  the  resulting  slipperiness.  Dry  surface 
conditions  improve  traction  and,  therefore,  little  attention  has 
been  paid  to  their  effect.  Nevertheless,  evaporation  from  moist 
soils  almost  invariably  results  in  a  surface  layer  that  is  sig¬ 
nificantly  drier  than  the  soil  at  depth.  In  order  to  assess  the 
effect  of  dry  and  wet  surface  conditions  on  mobility,  it  is  recom¬ 
mended  that  the  present  tire-soil  model  be  expanded  to  include 
layered  soil  conditions  as  well  as  soil  conditions  with  a  gradual 
increase  of  strength  with  depth.  Such  a  model  would  allow  the 
analysis  of  slipperiness  as  well  as  of  dry  surfaces  and  could  be 
used  to  develop  criteria  for  the  evaluation  of  field  tests. 

The  development  and  validation  of  tire-soil,  model  presented 
in  this  report  is  based  on  experiments  performed  in  either  purely 
fric'ional  or  purely  cohesive  soils.  Theoretically,  the  tire-soil 
model  does  not  differentiate  between  purely  frictional  and  purely 
cohesive  soils  since  the  soil  is  modeled  therein  by  its  strength 
properties  that  include  both  friction  and  cohesion.  Nevertheless, 
it  is  desirable  to  compare  performance  predictions  by  the  model 
with  experiments  and  make  modifications  in  the  model,  if  necessary. 
It  is  recommended  that  tire  performance  tests  in  a  varieLy  of 
c ohes ive- f r ic t iona 1  soils  be  conducted,  the  strength  properties 
of  these  soils  be  determined  and  the  data  made  available  for 
analysis  and  further  validation  of  the  model. 
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APPENDIX  A 


DETERMINATION  OF  SOIL  STRENGTH  PARAMETERS 
FOR  THE  PURPOSE  OF  MOBILITY  EVALUATION 


In  t  roduc  t ion 

Soil  strength  has  been  recognized  as  the  soil  property  that 
controls  the  behavior  of  soil-running  gear  systems  in  off-road 
locomotion.  Consequently,  in  the  tire-soil  model  presented  in 
this  report  the  soil  is  modeled  by  its  strength  properties.  For 
both  the  use  of  the  strength  properties  in  the  model  and  the 
characterization  of  soil  for  mobility  purposes  in  general  it  ! s 
essential  that  the  strength  properties  of  soil  and  the  conditions 
that  apply  to  the  laboratory  and  field  determination  of  these 
strength  properties  for  Lhe  purpose  of  mobility  evaluation  be 
do  f i ned . 

Definition  of  Soil  Strength 

Soil  strength  is  defined  as  the  maximum  shear  stress  that  can 
be  developed  on  a  failure  surface.  If  this  shear  strength  is  a 
linear  function  of  the  normal  stress  acting  on  the  failure  surface, 
then  the  following  equation  holds 

s  =  c  +  a  tan  .  (A-l) 

n 

where  c  and  c  are  the  Coulomb ian  shear  strength  parameters, 

lhe  linear  dependence  of  the  shear  strength  on  the  normal 
stress  is  an  approximation.  In  some  cases  it  may  be  necessary 
to  consider  nonlinear  strength  relationships;  the  use  of  non¬ 
linear  relationship  in  wheel-soil  interaction  problems  is  dis¬ 
cussed  in  Ref.  A.l.  For  the  purpose  of  the  following  discussions 
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it  is  assumed  that,  within  the  range  of  normal  stresses  of  interest, 
the  soil  strength  can  be  satisfactorily  approximated  by  the  lineai 
relationship  represented  by  Eq .  (A-l). 

The  Coulombian  soil  strength  parameters  are  not  material  char¬ 
acteristics  inasmuch  as  their  value  depends  on  Lhe  time  rate  of 
loading  and  the  drainage  of  the  soil  during  loading.  Also,  the 
Coulomb  parameters  refer  to  normal  and  shear  stresses  that  act  in 
a  plane  perpendicular  to  the  failure  surface  and  are,  therefore, 
essentially  two  dimensional  in  nature.  The  stress  and  strain  con¬ 
ditions  that  obtain  perpendicular  to  that  plane  may  affect  the 
value  of  the  Coulomb  parameters.  In  experiments  for  the  determina¬ 
tion  of  the  Coulomb  strength  parameters  for  mobility  purposes,  the 
transverse  stress  and  strain  conditions  should  duplicate  those  that 
exist  in  Lhe  field  in  the  tire-soil  interaction. 

The  following  discussion  is  intended  to  define  the  various 
conditions  in  tire-soil  interaction  that  have  a  bearing  on  the 
Coulomb  strength  parameters: 

o  Drainage  conditions  in  the  soil  are  generally  very  close 
to  the  hypothetical  undrained  condition  during  the  pas¬ 
sage  of  a  tire.  Exceptionally,  the  degree  of  drainage 
during  the  application  of  wheel  load  to  the  soil  may  be 
important  in  the  case  of  submerged,  relatively  well 
draining  soils,  such  as  fine  sand.  For  the  purpose  of 
general  mobility  evaluation,  undrained  conditions  may 
be  assumed  and  the  above  submerged  condition  treated 
as  an  exception. 

®  The  loading  rate  of  soil  in  tire-soil  interaction  is  di¬ 
rectly  proportional  to  the  translational  velocity  of  the 
tire.  Very  roughly,  the  soil  loading  rate  may  be  esti¬ 
mated  by  assuming  that  the  stresses  on  a  soil  element 
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r i s e  from  Lhcir  initial  stale  Lo  peak  during  a  Lime  in¬ 
terval  it  takes  the  tire  to  pass  half  of  the  length  of 
the  contact  area.  For  example,  a  Lire  traveling  at  the 
low  speed  of  i  ft/sec  and  having  a  1  ft  Jong  contact 
area  would  impart  fuJJ  loading  to  Lhe  soil  in  1/bth  of 
a  second.  Loading  rates  encountered  in  mobility  prob¬ 
lems  generally  fall  in  the  category  of  quick  or  rapid 
loading.  While  the  strength  of  quickly  loaded  soils 
may  be  considerably  different  from  that  obtained  with 
slow  loading  rates,  it  is  generally  not  too  sensitive 
to  changes  in  the  loading  rate  within  the  range  of 
translational  velocities  of  interest.  Thus,  in  strength 
testing,  loading  rates  that  are  approximately  in  the 
range  of  the  rate  anticipated  in  t he  field  are  acceptable 

•  Strain  conditions  in  tire-soil  interaction  are  generally 
Lhree  dimensional  since  the  compaction  of  Lhe  soil  be¬ 
neath  the  tire  is  generally  accompanied  by  lateral  dis¬ 
placement.  Thus,  even  though  a  two  dimensional  model  is 
used  for  the  computation  of  stresses,  Lhe  strength  paramo 
tors  of  the  soil  in  that  model  should  he  determined  by 
testing  methods  that  allow  for  strains  in  the  third 
dimension . 

The  development  of  soil  strength  is  always  accompanied  by 
volumetric  strain.  IC  the  displacement  constraints  in  a  particu¬ 
lar  problem  involving  failure  of  the  soil  arc  such  that  the  volu¬ 
metric  strain  associated  with  the  soil  strength  cannot  bo  attained, 
the  soil  strength  cannot  be  developed  and  the  use  of  a  so-called 
partially  mobilized  soil  strength  is  in  order.  In  soil  mechanics 
this  case  is  known  as  local  snear  failure,  as  contrasted  with  the 
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case  of  general  shear  failure  .hen  the  shear  strength  is  fully 
mobilized.  The  situation  in  tire-soil  interaction  problems  is 
different  from  that  encountered  in  conventional  bearing  capacity 
problems  in  that  with  the  advancement  of  the  tire  the  soil  is 
progressively  compacted  and  the  volumetric  strain  associated  wj Lh 
the  passage  of  the  tire  is,  in  Lbe  case  of  interest,  generally 
sufficient  to  mobilize  the  full  strength  of  soil.  Further  dis¬ 
cussion  of  strain  conditions  is  presenteu  in  connection  with  the 
use  of  triaxial  tests  for  the  determination  of  soil  strength 
properties . 


Laboratory  Determination  of  Soil  Strength  for  Mobility  Evaluation 

Two  types  of  testing  technique  are  used  widely  in  soil  mechanics 
practice  for  the  determination  of  soil  strength:  the  direct  shear 
type  and  triaxial  type  tests.  The  direcL  shear  tests  represent 
plain  strain  conditions,  whereas  in  the  triaxial  test  the  inter¬ 
mediate  principal  stress  equals  the  minor  principal  stress.  Strain 
conditions  in  the  triaxial  test  correspond  Lo  the  stress  conditions 
on  ■=  n~  .  This  latter  strain  condition  approximates  the  ones  in  the 

L  J 

field  more  closely  than  the  plain  strain  and,  therefore,  the  deter¬ 
mination  of  soil  strength  properties  by  triaxial  tests  is  more  ap¬ 
propriate  for  mobility  evaluation  than  direct  shear  tests.  The 
good  agreement  between  measured  interface  stresses  and  those  com¬ 
puted  on  the  basis  of  triaxial  shear  strength,  reported  in  Ref.  13, 
and  the  good  simulation  of  tire  performance  by  the  tire-soil  model 
using  strength  data  obtained  from  triaxial  tests  give  strong  sup¬ 
port  to  the  preference  of  triaxial  tests  for  strength  determination. 


The  samples  used  in 
obtained  in  the  field  or 


triaxial  tests  are  either  undisLurbed  ones 
ones  prepared  in  a  laboratory  mold.  In 
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principle,  undisturbed  samples  should  he  used  since  the  in  situ 
strength  of  soil  may  reflect  the  grain  structure  of  the  soil  that 
cannot  be  duplicated  in  laboratory  preparation.  This  is  particu¬ 
larly  true  for  residual  soils.  However,  undisturbed  sampling  of 
surface  soils  is  difficult  and  undesirable  where  field  variation 
in  soil  properties  requires  that  a  1 ?rge  number  of  samples  be 
taken.  Obviously,  field  testing  techniques  that  test  the  in  situ 
strength  of  soil  are  better  suited  for  the  purpose  of  mobility 
evaluation.  On  the  other  hand,  the  preparation  of  soil  beds  for 
tire  performance  tests  can  be  closely  duplicated  in  the  prepara¬ 
tion  of  soil  samples  for  triaxial  testing.  Thus  the  role  of  tri- 
axial  testing  in  mobility  evaluation  is  the  determination  of  strength 
properties  for  model  validation.  Triaxial  tests  are  also  necessary 
for  the  validation  of  field  testing  techniques. 

in  the  conventional  triaxial  testing  the  chamber  pressure 
(minor  principal  stress)  is  kept  constant  and  the  vertical  load  is 
increased  until  failure  occurs  or  (in  the  absence  of  a  clearly  de¬ 
fined  failure  point)  a  certain  strain  is  reached.  In  tire-soil 
interaction  the  minor  principal  stress  monotonieull y  increases  up 
to  1 L s  maximum  value,  then  decreases  monotonica 1 ly  throughout  the 
passage  of  the  wheel  —  as  can  be  seen  in  Fig.  A-l,  which  shows  the 
magnitude  and  direction  of  principal  stresses  in  the  failure  zones 
beneath  a  rigid  wheel.  The  stress  paths  corresponding  to  those 
principal  stresses  in  the  loading  and  unloading  phase  during  the 
passage  of  the  wheel  can  be  constructed  in  the  principal  stress 
space  (assuming  =  o^)  as  shown  in  Figs.  A-2a  and  A-2b  for  a 

rigid  wheel  and  in  Figs.  A-2c  and  A-2d  for  a  pneumatic  tire  0  In¬ 
dicated  by  light  lines  in  these  figures  are  the  equivolume  contour 
lines  that  can  also  be  obtained  in  triaxial  tests.  The  volumetric 
strains  associated  with  the  stress  path  of  loading  increase  together 
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Fig.  A- 1  Principal  Stresses  in  Soil  Beneath  a  Rigid  Wheel 

■with  the  increase  of  minor  principal  stress.  Thus  the  stress  and 
strain  conditions  that  duplicate  those  in  the  field  in  the  tri¬ 
axial  Lest  are  such  that  the  minor  principal  stress  (chamber  pres¬ 
sure)  wouid  have  to  increase  with  increasing  volumetric  strain. 
Since  conventional  triaxial  testing  techniques  do  not  allow  a 
a  variation  of  tire  chamber  pressure  during  loading,  appropriate 
modifications  were  made  to  the  triaxial  apparatus  available  at 
Grumman.  Description  of  the  testing  apparatus  ar.d  results  ob¬ 
tained  with  the  modified  technique  are  given  below. 

Triaxial  Testing  for  Mobility  Evaluation  Purposes 

In  order  to  apply  loading  rates  comparable  Lo  that  in  tire- 
soil  interaction,  the  triaxial  chamber  was  mounted  on  an  Instron 
TM4  loading  frame  that  allowed  a  strain  controlled  application  of 
the  vertical  load  at  rates  varying  from  0.05  to  5  cm/min.  Even 


the  highest  rate  oi  loading  thaL  tin*  loading  frame  allowed  was  less 
Lhan  that  corresponding  to  a  3  it/sec  travel  velocity;  the  re¬ 
sources  of  this  program,  unfortunately,  did  not  allow  for  a  major 
modification  of  the  loading  frame  .  Since,  in  mobility  problems, 
soils  are  generally  less  than  100  percent  saturated,  compression 
of  the  sample  occurs  during  loading,  and  Lhe  deformed  cross  section 
of  the  sample  cannot  be  determined  .ithout  measuring  the  lateral  de- 
formation  of  the  sample.  Initially,  one  circumferential  gauge  with 
a  centimct  'r  scale,  placed  at  mid-height  of  the  sample,  was  applied; 
it  was  found ,  ov’ever,  that  at  high  loading  rates  the  sample  deforma¬ 
tion  was  a  symmetrical  about  Lhe  mid-height,  tnd  it  was  necessary  to 
appiv  throe  c  i  rcu;  l  erential  otuges.  Since  i  t  was  noL  possible  to 
read  :  hesi  ,.au  tea  c  a  n  at  modera  t  •.  i  >adi  nr  rates ,  a  closed  circuit 
tel  ision  came  a  i  th  a  1  Fno  tape  recorder  was  used  to  record  the 
pi  sition  oi  ..  i  ga  . a  during  tin  test;  readings  were  made  by  re- 
playing  the1  video  t  -  and  stopping,  it  at  appropriate  intervals  to 
record  the  read .  in  picture  m  the  setup  is  shown  in  Fig.  A-3. 


a.  A-  :  arieule  Chamber  Pressure  Triaxial  Testing  with  Video 

:  pe  Re  ording  oi  tin  C  i  reumferer  Lia  1  Gauges 
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Triaxial  tests  with  the  above  technique  were  performed  on 
samples  of  Buckshot  clay  received  from  WES.  Figures  A-4  and  A- 5 
show  the  Mohr  circles  obtained  by  this  method  ror  tests  on  2.9x7 
in  samples  of  Buckshot  clay  prepared  at  14  percent  moisture  con¬ 
tent  and  Lested  at  a  loading  rate  of  5  and  0.5  cm/min,  re¬ 
spectively.  The  Mohr  circles  indicate  almost  identical  strength 
properties  for  these  two  loading  rates;  differences  between  the 
two  tests  are  within  the  range  of  accuracy  of  this  type  of  tcsL. 
Thus  at  the  loading  rates  applied  in  the  tests  strength  properties 
of  the  Buckshot  ciay  arc  not  significantly  affected  by  the  magni¬ 
tude  of  the  loading  rate. 

Unfortunately,  even  at  the  highest  lo.s  ,  i  ng  rate  '.hat  can  be 
employed  wiLh  the  loading  frame,  .bout  20  se  unds  elapse  before 
the  maximum  load  or  about  lO  percent  vertji  1  st’  lin  is  reached. 
Investigations  performed  with  a  dynami  loading  m.  bine  (Ref.  A.') 
showed  that  at  higher  loading  rates  the  effe  t  of  1  <  ling  nLt  be¬ 
comes  more  pronounced  and  a  gain  in  strength  wit  i  n  •  ase  of 
loading  rate  can  be  conclusively  observed. 

Mohr  circles  for  another  typical  test  on  KuekshoL  lay  are 
shown  in  Figv  A-b„  The  moisture  content  and  the  loauin 
the  sample  in  this  test  was  tie  same  as  in  test  Ill,  Inn  the  nasi  'in 
laLeral  pressure  of  about  0.5  kg/sq  cm  was  reached  a  5  percent 
strain  instead  of  at  10  percent,  as  in  test  Ill.  Tin  cohesion 
intercept  in  the  two  tests  are  about  the  same,  but  the  1  fiction 
angle  is  higher  in  test  B7  then  in  Bl. 

The  measurement  of  the  lateral  deformation  by  the  eir  uinler  n- 
tial  gauges  allows  the  evaluation  of  the  void  ratio  durin  i hi  test. 
Figure  A-7  shows  the  changes  in  void  ratio  as  a  function  1  lateral 
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Mo  hi'  Circles  for  Variable  Chamber  Pressure  Triaxial  Test 
B-l  -  Buckshot  Clry.  Moisture  Content  =  34%,,  Loading 
late:  5  cr  L 


0  1.5  2 

q  cm) 

Fig.  A- 5  Mohr  Circles  for  Variable  Chamber  Pressure  Triaxial  Test 
B-4  -  Buckshot  Clay.  Moisture  Content  =  34%,  Loading 
Rate:  0.5  cm/ min 
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Fig.  A-6  Mohr  Circles  for  Variable  Chamber  Pressure  Triaxial  Test 
BS-7  -  Buckshot  Clay.  Moisture  Content  =  34%,,  Loading 
Rate  =  0.5  cm/min 
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pressure  in  test  Bi.  A  relationship  between  void  ratio  and  stress 
states  in  the  slip  line  field  inav  be  established  by  these  measure¬ 
ments,  and  total  volume  change  associated  with  the  development  of 
the  slip  line  field  may  be  computed. 

Triaxial  tests  with  the  above  technique  were  also  performed 
on  materials  received  from  the  National  Tillage  Machinery  Labora¬ 
tory.  Figure  A-8  shows  the  Mohr  circles  for  the  Decatur  silty 
loam.  Moisture  contact  of  the  sample  was  8.7  percent  and  bulk 
density  1.2  g/cu  cm,  which  approximately  corresponded  to  that 
of  the  upper  2.5  in,  in  the  tire  performance  rest.  From  this 
test  the  Coulomb  parameters  are  c  =  130  lb/sq  ft  and  :  =  21°. 
Figure  A-9  shows  the  Mohr  circles  for  the  Vaiden  clay.  Moisture 
contact  of  the  sample  was  25  percent  and  the  bulk  density 
1.08  g/cu  cm,  approximately  corresponding  to  that  of  the  upper 
2.5  in.  in  the  tire  performance  tests.  The  Coulomb  parameters 
for  this  condition  were  found  to  be  c  =  300  lb/sq  ft,  T  =  12°. 

In  the  above  tests  an  al  tempt  was  made  to  duplicate  loading 
rates,  stress  paths,  and  strain  conditions  as  they  occur  in  the 
field.  Since  the  resources  of  this  program  did  not  allow  for 
major  modifications  of  the  equipment,  loading  rates  and  variation 
of  chamber  pressure  had  to  be  compromised.  It  is  believed,  never¬ 
theless,  that  even  with  the  restrictions  imposed  by  the  capabili¬ 
ties  of  the  apparatus,  the  test  results  are  acceptable  for  most 
practical  purposes.  An  advantage  of  the  proposed  testing  tech¬ 
nique  is  that  a  friction  angle  is  obtained  in  a  single  test,  whereas 
several  tests  are  required  for  the  determination  of  friction  angle 
with  the  conventional  method. 

For  further  development  of  the  testing  technique  additional 
information  on  stress-strain  conditions  in  tire-soil  interaction 
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Fig.  A  8  ‘l/1  hr  Circles  for  Variable  Chamber  Pressure  Triaxial  Test 

npn"-^Cat^r9Qilty  Loam-  Moisture  Content  =  8.7%,  Bulk 
uensi ty :  1.29  cu  cm 


Fig.  A  Mohr  Circles  for  Variable  Chamber  Pressure  Triaxia]  Test 
VI  -  V;  iden  Clay.  Moisture  Content  =  25%.  Bulk 

r>  i  f-»r  . 


is  required.  Such  information  will  be  forthcoming  in  connection 
with  the  research  project  aimed  at  the  determination  of  the  effect 
of  speed  on  tire  performance  where  velocity  fields  associated  with 
the  slip  line  fie]ds  will  be  determined  and  strain  fields  from  the 
velocity  field  may  be  derived. 

Concerning  the  testing  facility,  a  loading  frame  that  allows 
for  faster  loading  rates  is  necessary.  Chamber  pressure  regula¬ 
tion  that  would  allow  for  any  desired  variation  of  the  chamber 
pressure  with  vertical  strain  preferably  by  computerized  control 
is  another  necessity  for  the  development  of  a  triaxial  testing 
technique  geared  toward  the  needs  of  mobility  evaluation. 

Determination  of  Soil  Strength  Properties  by  Field  Tests 

Many  types  of  field  tests  have  been  designed,  and  several  of 
them  are  regularly  used  in  off-road  mobility  research  and  engineer¬ 
ing  for  the  characterization  of  soil  properties.  Some  of  these 
are  intended  to  measure  the  strength  properties  of  soils;  others 
are  used  to  obtain  some  measure  of  soil  strength  that  is  then  re¬ 
lated  to  running  gear  performance.  Every  field  instrument  that  is 
forced  into  'he  ground  and  records  a  displacement-force  relation¬ 
ship  yields  some  information  on  soil  properties;  i  ’  the  instrument 
displaces  the  soil  in  such  a  way  that  failure  zones  develop,  then, 
at  least  theoretically,  it  is  possible  to  evaluate  the  soil  strength 
parameters  from  the  results  of  the  field  test.  In  the  following  the 
suitability  of  present  field  techniques  for  the  determination  of 
soil-strength  properties  is  discussed,  and  methods  for  the  evalua¬ 
tion  of  these  field  tests  are  proposed. 
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The  Role  of  Volumetric  Strain  in  Field  Tests 


As  mentioned  before,  the  prerequisite  for  the  use  of  soil 
shear  strength  in  the  determination  of  stress  states  in  soil  by- 
plasticity  theory  is  Lhat  volume  changes  associated  with  the 
shear  strength  may  take  place  in  the  loading  process.  In  tire- 
soil  interaction  problems  of  interest  this  requirement  is  generally 
satisfied  since  the  volume  of  soil  displaced  by  the  tire  is  large 
compared  to  the  volume  of  soil  affected,  as  shown  in  Fig.  A-lOa. 

In  contrast,  in  a  cone  penetrometer  test  the  soil  volume  displaced 
by  the  cone  is  small  compared  to  that  volume  of  soil  affected  by 
the  failure  mechanism,  as  shewn  quantitatively  in  Fig.  A-lOb.  It 
should  be  noted  that  the  areas  shown  in  Fig.  A-lOb  are  cross  sec¬ 
tions  of  bodies  of  revolution;  to  obtain  volumetric  proportions 
the  areas  have  to  be  multiplied  by  the  respective  distances  of 
their  center  of  gravity  from  the  axis. 

Figure  A-lOc  shows  quantitatively  the  soil  volume  affected  by 
a  plate  sinkage  tost.  The  proportion  of  the  volume  of  soil  dis¬ 
placed  by  the  plate  to  that  affected  by  the  failure  mechanism  is 
greater  than  in  the  case  of  cone  penetration  and,  in  most  cases, 
is  sufficient  to  balance  the  volume  changes  associated  with  the 
development  of  failure  mechanism.  In  the  case  of  ring  shear  tests, 
the  volume  change  requirements  associated  with  the  failure  mechanism 
assumed  in  conventional  evaluation  procedures  are  satisfied,  but  in 
many  cases  that  failure  mechanism  is  not  valid.  A  more  detailed 
discussion  of  the  failure  mechanism  of  ring  shear  tests  is  given 
lacer  in  this  appendix. 
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b)  Field  Determination  of 
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Penetrometer 
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Fig.  A-IO  Comparison  of  Soil  Volumes  Displaced  and  Affected 
by  Failure  Mechanisms  in  s,b,  and  c 
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De terminal ion  of  Soil  Shear  Strength  Parameters  from  Cone 
Penetration  Tests 

In  the  tire  performance  tests  performed  at  WES,  cone  penetra¬ 
tion  tests  were  used  for  the  characterization  of  soil  properties. 

To  simulate  the  tire  performance  tests  by  the  proposed  tire-soil 
model  it  was  necessary  to  establish  the  Cou Lomb  soil  strength 
parameters  on  the  basis  of  cone  penetration  test  data.  To  this 
end  the  results  of  theoretical  analyses  as  well  as  triaxial  tests 
were  used. 

For  the  theoretical  analysis  of  cone  penetration  tests  a  com¬ 
puter  program  was  available  at  Grumman  that  computes  the  slip  line 
field  and  associated  stresses  as  well  as  the  volume  changes  in  the 
slip  line  field  for  various  cone  angles  assuming  shear  zones  up  to 
the  level  of  the  base  of  the  cone  and  a  surcharge  above  this  level 
for  the  consideration  of  depth  effect.  The  theoretical  basis  of 
this  program  and  other  details  are  described  in  Ref.  A. 3.  A  slip 
line  field  generated  by  thir'  program  assuming  strength  properties 
of  the  Buckshot  clay  is  shown  in  Fig.  A-ll.  Volume  change  computa¬ 
tions  using  void  ratio  [minor  principal  stress  relationships  deter¬ 
mined  by  triaxial  tests  (Fig.  A-7)  performed  for  the  case  when  the 
base  of  the  cone  is  level  with  the  surface]  showed  that  the  volume 
displaced  by  the  cone  was  generally  sufficient  to  produce  the  volume 
changes  necessary  for  the  lull  mobilization  of  soil  shear  strength. 
Thus  for  the  soil  condi ti<  ns  in  the  WES  tire  performance  tests  the 
cone  penetration  tests  may  he  regarded  as  a  measure  of  soil  strergth. 
It  should  be  noted,  however,  that  in  the  WES  tests  the  love'  t  rela¬ 
tive  density  of  the  Yuma  sand  was  about  50  percent  and  the  degree 
of  saturation  of  the  Buckshot  clay  about  95  percent.  At  that  rc  la- 
t  ive  density  of  sand  and  degree  of  saturation  of  clay  both  soils 
exhibit  relatively  small  volume  changes,  and  the  finding  thaL  for 
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the  soil  conditions  used  in  these  tests  cone  penetration  resistance 
is  representative  of  soil  strength  does  not  necessarily  mean  lhaL 
this  would  be  the  case  for  soils  exhibiting  larger  volume  changes. 

The  strength  properties  of  the  Yuma  sand  were  established  on 
the  basis  of  relative  density.  An  experimental  relationship  be¬ 
tween  the  relative  density  of  Yuma  sand  and  cone  index  gradient 
was  developed  at  WES  (Fig.  A-12).  Triaxial  tests  performed  at  WES 
on  Yuma  sand  were  also  related  to  its  relative  density  (Ref.  A. 4). 
On  the  basis  of  these  data  the  following  linear  relationship  be¬ 
tween  cot  :  and  relative  density  was  established 

cot  o  =  1.64  -  0.68  x  !;  .  (A-2) 

r 

This  relationship  is  valid  for  D  >  50  percent.  Figure  A-13  shows 
Eq.  (A-2)  and  the  values  of  :>  determined  by  various  triaxial 
tests.  The  friction  angles  represented  by  Eq.  (A-2)  were  close 
to  those  obtained  in  variablr  chamber  pressure  tests  and  by  theo¬ 
retical  analysis  of  the  cone  penetration  test  (Ref.  A. 3). 


For  the  Buckshot  clay,  empirical  relationships  between  the 
average  cone  index  and  strength  parameters  have  beet  developed  a! 
WES  (Ref.  A.  5).  These  relationships  have  been  compared  with  re¬ 
sults  of  variable  chamber  pressure  triaxial  tests  and  .heoretical 
analyses.  The- e  comparisons  showed  good  agreement  with  the  value 
of  cohesion  intercept  est 'blished  by  WES  but  indicated  a  somewhat 
higher  frictior  angle.  Accordingly,  the  following  relationship 
was  used  in  tire  performance  simulation: 


c (psi)  =  Cl/ 12 . 5 

T(degrees)  =  Cl/4 


C  A  -  3 ) 


Regarding  the  evaluation  of  soil  strength  from  cone  penetra¬ 
tion  test  data  the  following  comments  are  pertinent. 
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The  relationship  established  between  friction  angle  and  cone 
index  gradient  for  Yuma  sand  is  only  approximately  valid  Cor  other 
cohesionless  materials.  Theoretical  analyses  reported  in  Ref.  A. 3 
show  that  in  cohesionless  materials  the  cone  penetration  resistance 
depends  not  only  on  the  friction  angle  but  also  on  the  unit  weight 
of  material  and  the  interface  fricLion  Lhai  develops  on  the  surface 
of  the  cone.  For  a  particular  materiaL,  a r'  the  Yuma  sand,  there  is 
a  relationship  between  friction  angle  and  unit  weight  (expressed  by 
relative  density)  and,  therefore,  the  cone  index  gradient  is  a  mea¬ 
sure  of  strength.  For  another  cohesionless  material  this  relation¬ 
ship  may  be  different  and  the  cone  index  gradient  exhibits  a  dit- 
ferent  relationship  with  the  fricLion  angle  of  LhaL  material.  This 
theoretical  conclusion  is  supported  by  tire  performance  tests  per¬ 
formed  in  mor'ar  sand  (Rer.  L2 )  that  showed  a  pu 1  1  -coe f fir ien l  re¬ 
lationship  different  from  that  for  the  Yuma  sand,  although  the 
mortar  sand  is  in  many  respect.'-  very  similar  to  the  Yuma  sard.  Foi 
an  angular  sand,  or  a  sand  composed  of  grains  that  have  a  much  dll 
f event  specific  gravitv,  the  difference  would  be  even  ;rcater.  In 
addition,  in  very  loose  sane's  the  com  penetration  cannot  mobilise 
the  full  frictional  strength  and,  therefore,  it  is  not  a  unique 
measure  of  strength  in  loose  materials. 

The  relationship  established  between  Coulomb  strength  parame¬ 
ters  and  cone  index  for  the  Buckshot  clav  is  probably  an  acceptable 
approximation  foi  other  clays  that  are  saturated  close  to  100  per¬ 
cent  ,  but  may  not  be  valid  even  for  the  Buckshot  lay  at  a  dil  terenL 
degree  of  saturation.  Theoretical  analyses  show  that  the  unit 
weight  has  but  little  influence  on  thi  cone  penetration  resistance 
in  cohesive  soils  whose  friction  angle  is  very  low.  Alst  ,  there 
is  little  uncertainty  about  the  interface  adhesion  and  friction 
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Lhat  in  cohesive  soils  appears  to  be  fully  mobilized  in  cone  pene¬ 
tration  tests.  Thus  the  cone  index  is  a  lair  measure  of  soil 
strength  in  cohesive  soils  with  a  high  degree  of  saturation. 

Partially  saturated  cohesive  or  cohes ive - fric tional  soils  may 
exhibit  a  high  degree  of  compressibility.  In  such  soils  cone  pene¬ 
tration  cannot  mobilize  the  full  soil  strength  and,  therefore,  is 
not  a  good  measure  of  it.  A  relatively  Large  proportion  of  the 
one  resistance  is  due  to  the  adhesion  on  the  surface  of  the  cone 
and  the  effect  of  friction  angle  becomes  suppressed.  Thus  cone 
penetration  tests  are  not  best  suited  for  the  determination  of  the 
^tnngth  properties  of  such  soils. 

n  though  the  above  problems  raise  serious  questions  about 
the  suitability  of  cone  penetration  tests  f«.c  the  determination  of 
rrength  properties  of  frictional  cohesive  soils,  she  conve- 
of  field  testing  by  penetrometers  makes  their  use  for  that 
purpose  very  desirable.  It  is  evident  that  for  the  determination 
t.w  strength  parameters,  I  and  c,  at  least  two  parame- 
ol  the  field  test  are  needed  even  if  the  unit  weight  is  esti- 
ited  or  obtained  by  other  means.  The  penetration  resistance  aL 

■->  e  (C I  at  z  =  0)  and  the  rate  of  increase  of  this  re- 
(.granient)  would  yield  two  parameters,  but  the  insensi- 
30  apex  angle  cone  to  the  friction  angle  and  the 
i-  nt  displacement  to  mobilize  the  shear  strength  in  loose 
nder  these  two  parameters  unsuitable  for  strength  evalua¬ 
tion  tor-  ical  considerations  indicate  that  the  same  problems 

probab  would  not  arise  (or  would  occur  only  in  a  much  lesser 
degree)  with  90  or  120°  apex  angle  cones.  These  could  be  used 
either  in  on  junction  with  ihe  present  30°  cone  or  separately. 

A  theoretical  analysis  should  be  made  to  establish  criteria  for  the 
penetration  is  Lance  measurements  and  evaluation  methods  for  such 
co  nes  . 
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De  terming  Lion  of  Soil  S hear  Strength  Pa rameters  from  Plate  Sinkage 
Tests 

The  volume  of  soil  displaced  by  the  plate  in  plate  sinkage 
tests  is  gem  rally  sufficient  to  mobilize  the  fuJl  shear  strength 
of  soil  in  the  failure  zones.  Thus,  in  principle,  a  plate  sinkage 
tests  is  better  suited  for  the  evaluation  of  soil  strength  parame¬ 
ters  then  the  cone  penetration  test.  For  the  evaluation  of  soil 
strength  parameters  from  plate  sinkage  tests  it  is  necessary  to 
establish  the  point  at  which  failure  conditions  are  reached  (bearing 
capacity  failure)  and  beyond  which  the  pressure -sinkage  curve  repre¬ 
sents  plastic  equilibrium  at  the  particular  sinkage.  For  a  given 
soil  slip  line  fieLds  may  be  constructed  on  the  basis  of  plasticity 
theory  (Ref.  A. 6)  and  the  bearing  capacity  determined.  For  the 
evaluation  of  strength  parameters  from  the  results  of  plate  sinkage 
tests  the  computation  procedure  would  have  to  be  inverted;  since 
various  combinations  of  *.  and  c  could  yield  the  same  bearing 
capacity,  either  two  tests  with  different  plate  size  are  necessary 
or  a  relationship  between  the  rate  of  pressure  increase  beyond  bear¬ 
ing  capacity  and  friction  angle  .vould  have  to  be  established.  This 
latter  method  is  believed  to  be  feasible  on  the  basis  of  theoretical 
considerations.  In  summary,  plate  sinkage  tests  are  suitable  for 
the  evaluation  of  strength  properties  but  evaluation  procedures 
would  yet  have  to  be  developed. 

Deterrtr  iation  of  Soil  Shear  Strength  FaratneLe rs  by  Ring  Shear  Tests 

Ring  shear  tests  were  originally  designed  and  are  being  used 
for  the  determination  of  soil  strength  parameters  for  motility  pur¬ 
poses.  However,  the  evaluation  of  these  tests  are  objectionable 
because  the  failure  mechanism  in  these  tests  consists  of  cblique 
failure  surfaces,  as  discussed  in  Ref.  A. 7,  tha :  do  not  include 
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the  horizontal  surface  of  the  ring.  Thus  shear  stresses  referred 
to  this  horizontal  surface  are  not  representative  of  the  shear 
strength  of  the  soil  and  are  generally  less  than  that.  Unfor¬ 
tunately,  by  coincidence,  this  feature  of  the  ring  shear  test  re¬ 
sulted  in  lower  apparent  strength  and  yielded  more  realistic  trac¬ 
tion  values  than  the  true  strength  because  the  full  shear  strength 
at  the  running  gear-soil  interface  generally  cannot  be  mobilized. 

A  more  detailed  discussion  of  the  mobilization  of  shear  strength 
is  given  in  Ref.  13. 

For  the  ring  shear  test  to  yield  true  soil  shear  strength 
parameters  the  proper  failure  mechanism  would  have  to  be  considered 
in  the  evaluation  procedure.  Another  alternative  would  be  to  use  a 
sleeve  that  would  be  pushed  into  the  ground  around  the  ring  shear 
pl.  te  to  prevent  lateral  failure  along  the  oblique  failure  surfaces 
ano  force  the  failure  surface  to  develop  a*.  the  base  of  the  ring. 
Further  research  is  needed  to  develop  the  proper  technique  of  ring 
shear  tests  for  soil  strength  property  determination. 
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APPENDIX  B 


COMPUTER  PROGRAM  FOR  TEE  PREDICTION  OF  TIRE  PERFORMANCE 

Description  of  Computer  Program 

The  computer  prograi  listed  at  the  end  ol  this  appendix  per¬ 
forms  aii  computations  necessary  to  define  the  geometry  of  the  tire 
centerline  and  slip  line  fields  as  postulated  by  the  tire-soil 
model"  computes  the  interface  stresses,  and,  by  appropriate  inte¬ 
gration  of  these  stresses,  the  load,  torque,  ">nd  drawbar  pull 
values  for  initially  assumed  entry  and  rear  angles;  and  finds  the 
solution  for  a  given  load  bv  an  iteration  procedure  in  which  the 
entry  and  rear  angles  arc  appropriately  changed. 

The  computer  program  consists  of  two  putts:  the  main  program 
("KTIRE")  and  the  subroutine  ("s  Ff").  All  computation.1  necessary 
to  determine  the  cooidi nates  of  the  nodal  points  and  associated 
stresses  for  a  single  slip  line  ficic  are  performed  ia  the  sub¬ 
routine,  all  others  in  the  main  program. 

The  flow  diagram  for  the  main  prograi  is  shown  u  .  B-l. 

In  the  main  program  input  dat.  are  read  in  1'  om  a  fi.  ic  Input 
soil  properties  may  he  either  the  values  of  ,  am  .  or  cone 

penetration  data,  cone  Index  gradient  for  fri<  lional  soils, 
average  cone  index  far  cohesive  soils.  The  progi  rn  computes  the 
c,  : ,  and  .  values  corresponding  to  the  corn  penetration  data 
from  the  relationships  established  in  Appendix  A. 

From  the  input  data  the  program  computes  the  interface  fric¬ 
tion  angle  ,  the  separation  angle  j,  and  the  hypothetical 
normal  stress  q  ,  that,  is  the  normal  stress  lor  an  infinitely 
small  forward  field.  If  this  hypothetical  normal  stress  is  greater 
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than  the  limiting  stress  then  there  is  no  forward  slip  line 

field.  Depending  on  whether  there  is  a  forward  field  or  not  the 
flow  diagram  separates  in  two  branches  (see  Fig.  B-l). 

In  the  case  of  both  forward  and  rear  slip  line  fields,  the 
computation  starts  with  the  forward  field.  For  an  assumed  entry 
angle,  the  main  program  computes  the  coefficient  i  for  the 
logarithmic  spiral  and  calls  subroutine  SLFI  that  computes  the 
normal  stress  at  a  ^  pertaining  to  the  assumed  entry  angle. 

The  main  program  iterates  on  the  entry  angle  until  the  normal 
stress  q^  at  equals  p^  within  the  limit  of  tolerance  then 

calls  the  subroutine  for  the  computation  of  the  rear  slip  line 
field.  The  extent  of  the  rear  slip  line  field  is  determined  by 
the  condition  that  it  should  end  at  angle  where  che  inter¬ 
face  normal  stress  equals  p^.  Between  angles  ^  and  the 

interface  normal  stress  equals  p^  as  postulated  by  the  tire- 
soil  model.  Should,  however,  the  interface  normal  stress  compu¬ 
tations  in  the  rear  field  show  that  they  are  less  than  p^  —  even 
if  the  rear  field  extends  up  to  the  inner  end  of  the  forward  field 
( 'd)  —  then  p^  is  made  equal  to  q^  ,  the  normal  stress  at 
computed  from  the  rear  field,  and  the  computation  with  the  new  p^ 
value  is  repeated.  This  situation  may  occur  o'-  high  slip  values 
when  the  relatively  large  outward  directed  interface  friction  angle, 
6,  reduces  the  hearing  capacity  of  soil  and  the  resulting  normal 
stresses  are  lower  than  p^.  In  such  cases  tire  geometry  resembles 
to  that  of  a  rigid  wheel. 

In  the  last  part  of  the  main  program  load,  torque  and  drawbar 

pull  values  are  calculated  on  the  basis  of  the  interface  stresses 

obtained  from  the  two  slip  line  fields  for  the  assumed  .•  and  , 

r  d 

angles.  In  an  iteration  procedure  these  angles  are  changed  until 


B-3 


the  load  equals  the  design  load  within  the  limits  of  tolerance. 

In  the  case  where  there  is  only  a  rear  slip  line  field  the  flow 
diagram  shown  in  the  right  side  of  Fig.  B-l  holds. 

The  flow  diagram  for  the  subroutine  SLFI  is  shown  in  Fig.  B-2 . 
for  the  computation  of  a  forward  slip  line  field.  Designations  and 
grid  are  shown  in  Fig.  B-3.  The  rear  field  is  computed  in  a  simi¬ 
lar  way  except  for  the  condition  applied  to  the  determination  of 
the  inner  end  of  the  slip  line  held.  In  the  forward  field  the  slip 
line  field  ends  at  while  the  rear  field  ends  at  tae  angle 

where  the  normal  stress  equals  p^  . 

For  the  computation  of  the  coordinates,  principal  stresses  and 
their  directions,  the  plasticity  theory  for  soils  applies.  The 
differential  equations  of  plasticity  are  replaced  by  difference 
equations  and  are  solved  numerically.  The  solution  procedure  is 
described  in  detail  in  the  previous  report  on  rigid  wheel-soil 
interaction  (Ref.  1J)  and  elsewhere  (Refs.  B.l  through  B.3).  For 
brevity,  only  that  )art  of  the  solution  procedure  is  discussed 
here  that  is  different  from  the  rigid  wheel  problem,  the  geometric 
boundary  conditions  at  the  interface. 


For  the  rigid  wheel  the  geometric  boundary  conditions  at  the 
interface  requires  that  the  coordinate  points  x  and  z  must  lie 
on  a  circle  with  radius  R.  for  the  tire  these  points  must  lie  on 
a  logarithmic  spiral.  For  numerical  computations,  the  logarithmic 
spiral  is  approximated  by  a  polygon  that  is  tangential  to  it.  Log¬ 
arithmic  spirals  have  the  property  that  their  tangents  make  a  con¬ 
stant  angle  (r<  )  with  the  normal  to  the  radius  allowing  a  simple 

s  p 

computational  procedure  for  the  determination  of  the  polygon  side 
on  which  x  and  z  must  lie.  Accordingly,  the  following  equa¬ 
tions  apply  for  the  computations  of  stresses  and  coordinates  at  the 
boundary  of  the  active  zone. 
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of  Computations 


X  . 


X  .  ,  .  +  r<  Fx . 


+  '  (z,-  ~  Z.  .) 


i,j  1  +  )(jF  l-l, j  o  l-l, j-1  ov  i-i.j-1  i-l, J 


z  .  .  =  z .  .  .  +  F  (x .  .  ,  -x.  .) 

i,j  l-l, J-1  l-l, J-1  i,J 


(  B  —  1 ) 


o_.  ,  =■  o ,  ,  .  4-  2  •  tan  I  o_.  ,  ,  •  (1_  ,  -  •  ,  .)  +  yC 


i,j  i-l,j 


i~l , J  v  i,j  1-1)1 


where 


o 


=  cotan  ) 

v  1>J"1  i ;  J 


C  =  z.  .  -  z.  ,  .  -  tan  9  •  (z.  .  -  x.  ,  .) 

i , J  i-l,  J  i,3  1-1,1 


F  =  tar.  i.  ,  .  i  -  i 

i  - 1 ,  1  - 1  s  p 


There  are  limitations  to  the  curvature  of  the  logarithmic 
spirals  that  are  incorporated  in  the  program.  These  limitations 
are 

In  the  front  field,  the  1  coefficient  is  limited 
to  such  a  value  that  the  z  coordinates  of  the  spiral 
be  always  positive  (below  the  surface). 

In  the  rear  field,  the  angle  of  inclination  to  the 
horizontal  of  the  side  of  the  polygon  at  a  particular 
point  is  limited  to  the  angle  of  a  straight  line  from 
that  point  to  the  end  point  of  the  forward  field  (no 
double  curvature  in  the  center  portion  is  allowed)  . 

If  both  limitations  hold  the  cent-v  line  geometry  reverts  to  ^hat 
of  a  tire  on  a  rigid  surface. 

The  flow  diagrams  in  Figs.  B-l  and  B-2  show  the  genera L  logic 
of  the  computations.  There  are,  however ,  some  provisions  in  the 
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programs  char  for  clarity  have  not  been  included  in  the  flow  dia¬ 
gram.  These  are  discussed  below. 

The  finite  difference  approximations  of  the  differential  equa¬ 
tions  of  plasticity  become  indeterminate  if  either  q>  or  c  equals 
zero.  Instead  of  providing  an  alternate  and  complicated  algorithm 
for  these  specific  cases,  a  minimum  value  of  1°  is  assigned  for 
T  and  0.2  lb/sq  ft  for  c.  The  effect  of  these  assignments  on 
performance  calculations  is  minimal. 

In  Ref.  13  it  was  shown  that  the  maximum  value  of  the  angle  of 
interfere  friction  is 

6  =  arc  tan  (sin  cp)  .  (B-2) 

At  this  value  of  the  interface  friction  angle  some  computa¬ 
tional  procedures  become  inaccurate  or  indefinite.  To  avoid  this 

computational  problem  a  maximum  value  of  0.99  6  has  been 

max 

assigned  to  6. 

In  the  various  iteration  procedures  that  are  employed  to  meet 
certain  equality  conditions  tolerance  limits  have  been  assigned  to 
each  condition,  on  the  basis  of  experience  gained  with  running  th. 
program  for  a  variety  of  conditions.  These  tolerance  limits  are 
consistent  with  the  degree  of  accuracy  desired  in  tire  performance 
predictions.  There  are  inevitable  inaccuracies  in  prediction 
theories  due  to  the  various  assumptions,  inaccuracies  in  the  de¬ 
termination  of  soil  properties  anti  the  approximations  inherent  in 
numerical  methods.  Thus,  generally  there  is  no  point  in  setting 
lower  tolerance  limits  in  the  equality  conditions  and  performing 
the  ensuing  larger  number  of  iterations  to  achi?ve  an  apparent 
accuracy  in  the  prediction  when  inaccuracies  from  other  sources 
would  dominate. 
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Convergence  of  Lhe  various  iterations  used  in  the  computer 
program  has  been  generally  established  in  the  numerous  runs  that 
were  required  to  develop  the  program  and  compare  predictions  with 
test  results.  Nevertheless,  it  is  conceivable,  that  for  some  un¬ 
foreseen  combination  of  variables  an  iteration  would  not  satis¬ 
factorily  converge.  A  time  limit  of  >00  •  conds  has  been  set 

for  the  program;  should  a  solution  be  not  reached  within  this  Lime 
limit  the  originator  of  the  program  should  be  consulted. 


Guidance  for  Users  of  the  Computer  Program 

Complete  listings  of  the  computer  programs  are  given  at  t  lie 
end  of  this  appendix. 

Designation  of  variable  names,  instructions  for  the  prepara¬ 
tion  of  data  files  and  other  comments  useful  for  'unning  the  pro- 
g  r  a  m  •  i  r  e  g  i  v  c  n  b  e  1  o  w  . 

Input  data  are  read  from  input  data  file  "TYR.DAT"  in  the  fol¬ 
lowing  order  and  2F  format. 


Com  Index  gradient  (pci)  LCGR 

Nr.iiinal  tire  radius  (ft)  _  RA 

Tire  load  (lb)  L <t)  ' 

Deflection  coefficient  DE  ] 

Slip  parameter  (j  )  [ S J 

Cohesion  in  forward  field 
(Ib/sq  ft)  [CF] 

Friction  angle  in  forward  field 
(degrees)  FF 


Cone  index  (psi)  iCIj 
Tire  width  (ft)  j 
Limit  pressure  (p,i)  PfD 
Slip  [SL 

Slip  parameter  K  [  SKF. 

Cohesion  in  rear  field  [ CR ] 

Friction  angle  in  rear  field  FR ] 


Unit  weight  in  forward  field 
(lb/cu  ft)  CF  i 


Unit  weight  in  rear  field  1 GR ] 


1 1  w  'j.  ,  '■!  VTPV“-  "TM1 


.'ju^nr 


Des Ignat  tors  of  the  input  variable'  in  thi  omputei  program 
are1  indicated  in  brackets. 

Cohesion,  friction  angle,  and  unit  weight  are  read  in  a  input 
variables  only  if  both  CGR  =  0  and  Cl  =  0.  If  t  hi  progran  is  L» 
be  used  with  cone  penetration  daLa,  the  value  of  ilher  CGR  or 
Cl  should  be  entered  in  the  data  file,  the  other  value  must  be  zero. 

Subscripted  Variables  and  Dimension  Statements 
Main  program  KTIRE 

J  designates  a  location  at  the  interface 

HH  ( J)  =  ... 

.1 

QQ(J)  =  q- 

= 

UU(J) 

z. 

J 

Triple  letter  subscripted  variabb  are  in  ior  sua  age  of  double 
letter  arrays. 


Subrou  t  i.ne 

SLI'I 

X  ( 1 ,  J  )  =  x 

i  J 

X(I,J  -  z 

i » J 

S(I,J)  =  ■ 

i  » j 

T  C I  >  J  )  =  ■ 

i  ,  j 
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The  dimension  statement  corresponds  to  a  48  x  16  grid 
1  -  1  to  49,  J  -  1  Lo  17)  for  the  computation  of  the  geome¬ 
try  of  slip  line  fields.  All  J  locations  are  not  necessarily 
used  in  actual  computations,  as  indicated  before.  Several  hundred 
computations  were  performed  satisfactorily  with  the  anove  grid 
size  and,  there  fore,  no  need  for  the  change  of  the  grid  size  is 
ant icipa  ted . 

Variables  with  one  subscript  are  as  follows: 

D(J)  =  h. 

H(J)  = 
q(.J)  -  q. 
n(.J)  =  ; 

!'(J)  =  x. 

(J)  =  z . 

1 

A ( .  1  ,  bid),  C(J)  are  auxiliary  variables. 

The  dimensions  of  the  variables  with  one  subscript  are  tied 
to  the  J  dimensions  of  the  variables  with  tw.i  subscripts. 

Limit 

In  t hv.  main  program  the  following  limits  are  set  for  the  entrv 


real 

■  angles . 

e 

(ARE) 

max 

=  70 

t 

r 

(ARMAX) 

ma  x 

ii 

■O 

O 

3 

r 

.  (ARMIN) 
mm 

=  0° 

B-L  l 


'  r?  *"  ■',‘lP7*r?7^rr?  • 


A  limit  of  20  was  set  for  the  numbe1*  of  iterations  (ITER) 
allowed  to  find  the  solution  for  a  given  load. 

A  limit  for  the  tire  radius  at  the  rear  entry  angle  (RAMIN) 
was  set  to  correct  tire  geometries  that  would  have  resulted  in 
negative  sinkages. 

Other  Constants 


The  initial  rear  angle  (AR)  is  set  at  10°. 

The  angle  (AD)  is  set  at  5°. 

A  constant  DF  defines  a  linear  variation  oi  6  over  the 
interface.  In  the  present  program  &  is  assumed  to  be  uniform 
over  the  interface.  For  this  case  a  value  of  DF  =1  is  assigned. 

In  the  main  program  XXX,  in  the  subroutine  XX,  is  used  to  define 

che  conditions  for  which  the  slip  line  field  is  computed  according 
to  the  following  schedule 

XXX  =  -  1  for  forward  slip  line  field 

XXX  =  +  1  for  rear  slip  line  field 

XXX  =  +  2  for  rear  slip  line  field  only 


Ou  t  pu  t 

The  following  results  of  the  computation  are  printed  out. 

Pull  coefficient  (PUN) 

Torque  coefficient  (TON) 

Sinkage  (SNK)  (inches) 

Note:  The  tire-soil  model  has  been  developed  for  driven  tires 

where  the  applied  torque  determines  the  direction  of  shear  stresses 


that  is  the  same  along  the  whole  interface.  In  the  case  of  towed 
tires  the  direction  of  shear  stresses  changes  over  the  interface 
so  that  the  resultant  torque  from  the  shear  stresses  is  zero. 
Therefore,  the  tire-soil  model,  as  presented,  should  not  be  used 
for  the  computation  of  the  drag  force  of  towed  tires. 


Reier  jnces 

B.l  Karafiath,  L. ,  "Plasticity  Theory  and  Stress  Distribution 

beneath  Wheels,"  Journal  of  Terramechanics ,  Vol.  8,  No.  2,  1971. 

B.2  Karafiath,  L.  and  Nowatzki,  E.,  "Stability  of  Slopes  Loaded 
over  a  Finite  Area,"  Highway  Research  Board  Record  No,  323, 
November  1970, 

B.3  Karafiath,  L. ,  "On  the  Effect  of  Pore  Pressures  on  Soil-Wheel 
Interaction,"  Proceedings  Fourth  International  Conference  for 
Terrain  Vehicle  Systems,  Stockholm,  April  1972. 
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COMPUTER  PROGRAM  LISTING 
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DIt’P.F  I  P’J  yfAF,  17  ),7(AC,  !7),F(49,17),T(AQ,  17), 
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IMPLICIT  17  EAL  (L)  v 

L  EL  cr  9,  F9  )=AS  IN  CS  110  CD9  )/F  ItJ  (  F9  )  )  /A 

r  IF  A  CL  9,  F  G  )  =  C  CS  (C9)+F0-T  (CCF  (F9)**2-CLS  ( F  c ) *  *2  ) 
PUP ( C9,  F9 ) =C CF ( C9 ) -SOFT (CCS  CD 9 )**2-CCF  ( F9 ) *  * 2  ) 
FPP(FG,TG,TP  )  =  FX7"  (  2  *  ( T  9  -  T  P  )  *  fi  I M  (F9)/CCF  (FG  )  ) 
TAD  (  FG  )-c  IM  (  FG  )  /C  Cc  (  F9  ) 

F  IF  ( A 9 j  F 9 ,  T9, TP,  F9  )=F0*  (A G  D9)*EXr  (  (T9-TP  )x 
+  F-  ID  (  ;  °  )/C  OF  (  F9  )  ) 
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